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INTRODUCTION 


The Institute for Computer Applications in Science and Engineering (ICASE) is oper- 
ated at the Langley Research Center (LaRC) of NASA by the Universities Space Research 
Association (USRA) under a contract with the Center. USRA is a nonprofit consortium of 
major U. S. colleges and universities. 

The Institute conducts unclassified basic research in applied mathematics, numerical 
analysis, and computer science in order to extend and improve problem-solving capabilities 
in science and engineering, particularly in aeronautics and space. 

ICASE has a small permanent staff. Research is conducted primarily by visiting scientists 
from universities and from industry, who have resident appointments for limited periods of 
time, and by consultants. Members of NASA’s research staff also may be residents at ICASE 
for limited periods. 

The major categories of the current ICASE research program are: 

• Numerical methods, with particular emphasis on the development and analysis of basic 
numerical algorithms; 

• Control and parameter identification problems, with emphasis on effective numerical 
methods; 

• Computational problems in engineering and the physical sciences, particularly fluid 
dynamics, acoustics, and structural analysis; 

• Computer systems and software for parallel computers. 

ICASE reports are considered to be primarily preprints of manuscripts that have been 
submitted to appropriate research journals or that are to appear in conference proceedings. 
A list of these reports for the period April 1, 1991 through September 30, 1991 is given in the 
Reports and Abstracts section which follows a brief description of the research in progress. 


^Presently, ICASE is operated at NASA Langley Research Center, Hampton, VA, under the National 
Aeronautics and Space Administration, NASA Contract No. NASl-18605. In the past, support has been pro- 
vided by NASA Contract Nos. NASl-18107, NASl-17070, NASl-17130, NASl-15810, NASl-16394, NASl- 
14101, and NASl-14472. 
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RESEARCH IN PROGRESS 


Saul Abarbanel 

The study of the stability of hyperbolic initial-boundary value problems has led to an 
ICASE Report by Carpenter, Gottlieb and Abarbanel entitled: "The Stability of Numerical 
Boundary Treatments for Compact High Order Finite Difference Schemes”. 

The above study relied heavily on numerical evaluations of various stability criteria. To- 
gether with David Gottlieb, additional theoretical studies were carried out to discover under 
which norms one might be able to prove stability. For 1st and 2nd order numerical boundary 
conditions this issue has been resolved and stability proved. This work now continues with 
regard to higher order boundary conditions and for the more general case of systems. 

The work (together with doctoral student J. Danowitz) on non-reflecting boundary con- 
ditions for compressible viscous flow past two-dimensional finite bodies continues. 


Remi Abgrall 

We have tried to design an Essentially Non Oscillatory reconstruction for functions de- 
fined on finite-element type meshes. Two related problems are studied : the interpolation 
of possibly unsmooth multivariate functions on arbitrary meshes and the reconstruction of 
a function from its average in control volumes surrounding the nodes of the mesh. Con- 
cerning the first problem, we have studied the behavior of the highest coefficients of the 
Lagrange interpolation functions which may admit discontinuities of locally regular curves. 
This enables us to choose the best stencil for its interpolation. The choice of the smallest 
possible number of stencils is addressed. Concerning the reconstruction problem, because of 
the very nature of the mesh, the only method that may work is the so called reconstruction 
via deconvolution method. Unfortunately, it is well suited only for regular meshes as we 
show, but we also show how to overcome this difficulty. The global method has the expected 
order of accuracy but is conservative up to a high order quadrature formula only. 

Some numerical examples which demonstrate the efficiency of the method are given in a 
forthcoming ICASE report. 


H. T. Banks and Fumio Kojima 

Work is continuing on the study of the corrosion detection problem arising in thermal to- 
mography. In conjunction with W.P. Winfree (Nondestructive Evaluation Sciences Branch, 
LaRC), we have developed practical corrosion detection algorithms based on nonlinear least 
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square estimation. Computational efficiency is a very important factor in implementing 
nondestructive evaluation techniques. The latest version of our algorithm is at least four 
times faster than the earlier one. The choice of the regularization parameter in our geomet- 
rical heat inverse problem has been studied also. A graphical method to find the optimal 
choice of the regularized parameters has been investigated; this method involves plotting the 
seminorm of the solution versus the residual norm. 


H. T. Banks and Ralph Smith 

Work is continuing on the development of a computationally viable infinite-dimensional 
control and identification methodology for the stable and optimal design of noise suppression 
systems. The basic approach is to combine approximation theory with time domain state 
space modeling to develop convergent computational algorithms for LQR control designs. 

The model to be considered consists of an exterior noise source which is separated from 
an interior chamber by an active wall or plate. This plate transmits noise or vibrations from 
the exterior field to the interior cavity via fluid/structure interactions thus leading to the 
formulation of a system of partial differential equations consisting of a 3-D acoustic wave 
equation coupled with elasticity equations for the plate. The control is implemented in the 
model via piezoceramic patches on the plate which are excited in a manner so as to produce 
pure bending moments. The validity of the model is being tested both numerically and with 
data from experiments by H.C. Lester and R.J. Silcox (Acoustics Division, LaRC). 

As an initial test of the control methodology, a simplified but typical model consisting 
of a 2-D interior cavity with an active beam at one end has been considered. We have 
found that a Legendre basis in the cavity coupled with a cubic spline discretization of the 
beam provides a scheme which efficiently handles the numerical difficulties resulting from 
the lack of a uniform margin of stability and the fact that the system is only weakly damped. 
Numerical tests have shown that when the feedback controls are implemented, substantial 
reductions in acoustic pressure can be realized with physically realistic voltages into the 
patches. Further numerical tests are being made to determine optimal patch placements, 
numbers and lengths for the 2-D model. 

Current and future work is also being directed toward developing an accurate, efficient 
and stable means of discretizing the 3-D model in which the interior cavity is taken to be 
cylinder with a circular active plate and sectorial patches. Once the 3-D approximation 
scheme has been completed, the identification of physical parameters and control of acoustic 
pressure will be implemented using techniques similar to those developed for the 2-D problem. 
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Alvin Bayliss 

In collaboration with L. Maestrello (Fluid Mechanics Division, LaRC) and A. Frendi 
(Vigyan, Inc.), we consider the evolution of unsteady disturbances in a supersonic boundary 
layer over a flexible surface. The problem is solved numerically. The unsteady flow field 
disturbances are computed by direct solution of the Navier-Stokes equations. The motion 
of the wall is computed by solution of the equations governing a flexible surface clamped 
at each end. These equations are coupled, as the plate equation is forced by the pressure 
computed from the Navier-Stokes equations which in turn requires the solution of the plate 
equation as a boundary condition. 

We consider the case where the plate is forced by an instability wave in the flow field 
as well as the case where the flexible source is excited by an acoustic disturbance. Results 
obtained to date demonstrate that the growth of disturbances can be enhanced when the 
bounding surface is flexible. 


H. Scott Berryman 

In cooperation with Intel Scientific Supercomputing Division, and with the collaboration 
of W. Gropp (Argonne National Laboratory) and R. Littlefield (Pacific Northwest Labora- 
tory) the behavior of the Touchstone DELTA system interprocessor communication subsys- 
tem has been studied. The goal is to gather data to aid in protocol development, application 
optimization and the development of interprocessor global operations such as the global 
reduction, and global transpose. 

In collaboration with S. Bokhari algorithms for the global transpose for two dimensional 
grids in general, and the DELTA system in particular, are being studied. The global trans- 
pose is needed in such scientific codes as multidimensional FFTs, unstructured mesh codes 
(for address translation using lookup tables) and ADI methods. 


Nicholas Blackaby 

Work on the stability of hypersonic flow over a flat plate in collaboration with S. Cowley 
(Cambridge University) and P. Hall was completed. A paper based on an earlier ICASE 
report by the authors has been submitted to the Journal of Fluid Mechanics. 

Research into vortex-wave interactions in compressible boundary layer flows has con- 
tinued. The weakly non-linear interaction of Tollmien-Schlichting waves and longitudinal 
vortices has been considered, in collaboration with F.T. Smith (UC London). The inviscid 
stability of a small-wavelength, fully nonlinear Gortler vortex state in a large-Mach number 
Blasius boundary layer was investigated during the LaRC/ICASE Transition and Turbulence 
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workshop; a report on this work has been submitted to the workshop proceedings. 

Together with M. Choudhari (High Tech. Inc.), the nonlinear development of longitudinal 
vortices in the presence of wall-curvature and cross-flow is presently being considered. In 
most practical aerodynamic applications there will be a significant amount of cross-flow 
normal to the principal direction of flow. Cross-flow has previously been shown to have a 
dramatic effect on the linear stability properties of Gortler vortices. We have found that, 
above a certain size of cross-flow, neutrally stable longitudinal vortices are possible over a 
convex surface. We are currently investigating the weakly nonlinear development of these 
modes (for both concave and convex wall curvatures). The problem is mathematically very 
similar to studies into the stability of stratified flows in geophysical applications. Thus we 
are able to apply some of the methods of the latter to our study. A Landau amplitude 
equation will be derived from which nonlinear development properties can be deduced. An 
ICASE report describing this work should follow early in 1992. 


Percy Bobbitt 

A paper concerning the background of, and results from, the hybrid laminar flow control 
(HLFC) experiment carried out in the LaRC’s 8-ft Transonic Preserve Tunnel on a swept, 7- 
foot chord model has been written and is awaiting final review and publication. A shortened 
version of this paper has also been prepared for presentation at the Annual AIAA Conference 
in Reno, Nevada and is now being reviewed. Both papers discuss the design of the HLFC 
airfoil’s upper-surface contour and the associated pressure distribution. Results for variations 
in Mach number, Reynolds number, level and extent of suction are illustrated and analyzed. 
The results are given primarily as plots of transition location, suction mass-flow requirements, 
pressure distribution and drag. 


Shahid Bokhari 

A major overhead in interprocessor communications on distributed memory machines is 
the time required to permute data blocks when using a ’combining’ algorithm. The time 
required to shuffle blocks in the ’standard exchange’ algorithm for complete exchange on hy- 
percubes is an important (though not the only) example. I am developing an enhanced DMA 
unit that will eliminate the overhead of this data permutation at the cost of a modest invest- 
ment in hardware. Since this overhead can be as much as 50% of the total communication 
time, this enhancement is likely to be of considerable practical value. 
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A theoretical analysis of the multiphase algorithm for complete exchanges on the hy- 
percube is also underway. In this case I am attempting to develop a method for efficiently 
obtaining the optimal algorithm, given a specific block size and hypercube dimension. 


Shahid Bokhari and Scott Berryman 

Work is underway on developing a class of algorithms for complete exchange on the 
Touchstone Delta Mesh. The complete exchange is at the heart of many important applica- 
tions. The mesh interconnect is much sparser than the hypercube and there is no possibility 
of an asymptotically optimal ’direct’ algorithm, as was the case for hypercubes. We are 
developing a multiphase combining algorithm which trades off the startups of the direct 
algorithm with the overhead of data permutation. A theoretical analysis is nearly complete 
and we are in the process of obtaining actual run times on the Touchstone Delta machine. 


Kurt Bryan 

Work is continuing with W.P. Winfree (Nondestructive Evaluative Sciences Branch, 
LaRC) on a thermo-mechanical inverse problem, specifically, using the surface mechani- 
cal response of a structure to external heating to determine its internal properties and/or 
locate defects. The focus is on structures which are thermally and mechanically isotropic 
and homogeneous. Analysis of the forward problem leads to a system of partial differen- 
tial equations, the heat equations coupled to the equations of linear elasticity. Due to the 
complicated nature of the equations, the computational solution of the inverse problem will 
probably rely on an optimization/fit-to-data approach. This requires an economical scheme 
for solving the forward problem. A program has been written to solve the forward problem, 
based on the method of boundary integral equations. This approach should be particularly 
efficient, since the inverse problem requires the mechanical response of the structure only on 
its surface. We are now pursuing appropriate convergence results for the numerical method 
and a sensitivity analysis to determine the optimal experimental configuration for locating 
internal defects. 

Work is also underway with M. Vogelius (Rutgers University) on a version of the inverse 
conductivity problem. The goal is to recover information about multiple internal cracks in 
an electrical conductor based on voltage and current flux measurements on its boundary. 
Theoretical results concerning the detection and identification of multiple internal cracks in 
a conductor have been proved. A numerical algorithm is also under development, based on 
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an integral equation formulation of the problem. This formulation and code should lead ] 

to results for the optimal electrode arrangement and computational solution of the inverse I 

problem. 

Thomas W. Crockett 

The iPSC/860 version of the parallel renderer was modified to generate sequences of 
images using both temporal and spatial compression. This is an attractive technique for 
producing live animation on the iPSC/860 given the low bandwidth of its external connec- 
tions. A companion interpreter for the compressed image format was developed which uses 
the X Window System and Unix sockets to display the output remotely on users’ work- 
stations. The interpreter will also run directly on the iPSC/860 to play back previously 
recorded sequences stored on the Concurrent File System. The combination of the renderer 
and interpreter serves as the first demonstration at Langley of the iPSC/860’s backend CIO 
Ethernet interface. The techniques employed have wide applicability in other distributed 
codes which are being developed for the machine. 

A new release of the iPSC/860 system software from Intel was installed. The new soft- 
ware, while providing critically needed functionality, proved to be less stable than the pre- 
vious version, resulting in a series of system crashes with major filesystem damage. Con- 
sequently, disk backup and restore procedures were refined so that recovery time could be 
reduced from 14 hours to 5 hours. Pressure from ICASE and other sites has convinced Intel 
to provide monthly software updates. 

Intel’s system software currently provides little assistance in sharing resources fairly 
among users. With our relatively large user community, this was becoming a serious problem. 

After investigating several alternatives, we were able to acquire a job scheduling package un- 
der development at the San Diego Supercomputer Center. We are currently acting as a beta 
site for this software, which we expect to be incorporated into a future release from Intel. 

The SDSC software has been very successful, virtually eliminating the resource management 
problem. 

Wai-Sun Don 

In collaboration with A. Solomonoff, the research about the accuracy and speed in com- 
puting the Chebyshev collocation derivative was completed, and a paper describing the 
results was submitted to the SIAM Journal on Scientific and Statistical Computing for pub- 
lication. 

Collaborating with P. Drummond (Fluid Mechanics Division, LaRC), a three dimensional 
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spectral code for solving the hydrogen vortex jet interaction with an oblique shock using 
Chebyshev collocation methods is currently under development. The preliminary results of 
the air vortex jet interaction with an oblique shock look quite promising. Chemical reactions 
with different species will also be included in the calculation to simulate the physics in the 
scramjet engine. 


Peter W. Duck 

The interaction between a shock wave and various forms of freestream disturbance (i.e. 
vorticity waves, entropy waves and acoustic waves) is currently under investigation. On 
passing through the shock, any one of these waves will, in general, generate all three classes 
of wave. The specific configuration under study involves a wedge (which is responsible for 
the creation of the shock), and the effect of disturbances ahead of the shock on the boundary 
layer on the wedge surface is of particular interest. There axe a number of other issues related 
to this problem currently under study, including the effect of wedge surface distortion, and 
the possibility of altering the shock strength. 


Thomas Eidson 

An explicit, finite-difference code developed by G. Erlebacher to study 3-D, compressible 
turbulence was converted to the Intel iPSC/860 to evaluate various performance program- 
ming ideas. Two versions which have different distribution strategies have been implemented. 
One version uses across-processor transposes to get the appropriate data on each processor 
to perform the tridiagonal matrix solutions. The second version solves the tridiagonal ma- 
trix problem with an algorithm that works with the appropriate data distributed across 
several processors. Both versions attain between 5 and 8 Mflops per node. Problem sizes 
up to 64x64x64 have been tested. Larger sizes will be tested when the initial condition code 
is running on the hypercube. More extensive data analysis, which includes an FFT algo- 
rithm, is also being implemented. This code has also been used to study several problems 
in programming in a multi-machine environment. The code was modified so that minimum 
differences exist between the CRAY version and the hypercube version. The code was also 
modified to run different segments of the code on different computers: the control segment 
runs on a workstation, the main calculations run on the hypercube and the monitor/graphics 
segments runs on a graphics workstation. 

A vector library (written in assembler) was written and tested. There does not appear 
to be any significant benefit of using a vector library over the compiler in the general case. 
This is partially due to significant improvements in the compiler. For code sections for which 
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the data could reside in cache and be reused, assembler coding may provide some benefit. 
An example would be an FFT. 

A project to develop a code transformation tool was begun. With an increasing number 
of computers with varying architectures becoming available, moving codes between comput- 
ers is being done more frequently even though it is becoming more difficult. The design of 
new compilers and automatic pre-processors is not keeping up with the demand. The trans- 
formation tool is designed to aid the user in porting codes by handling the details of the code 
changes with fewer errors, thus reducing the debugging effort. The user will be required to 
determine the appropriate transformations and thus do the appropriate data dependency or 
data flow analysis. This strategy will hopefully provide a useful porting tool more quickly 
than developing a compiler or pre-processor that does most of the code analysis. 


Gordon Erlebacher 

Current collaboration with S. Sarkar is focused on the analysis of compressible shear 
turbulent flow databases that he generated in the last year. The objective is to elucidate the 
effect of compressibility on the well-known structures of incompressible shear-flow turbulence. 
The results of this work will be submitted to Physics of Fluids and presented at the APS 
meeting. 

First hand knowledge of wavelet theory and its potential applications to CFD was ob- 
tained at a meeting on wavelets at Princeton University. This led to a collaboration with 
C. Streett (Fluid Mechanics Division, LaRC), in which several diagnostics were developed 
based wavelet transforms to perform a quantitative study of long time signals obtained from 
the direct simulation of finite-length incompressible Taylor Couette flow. This allowed the 
frequency content of the signal to be analyzed as a function of time. The results of this work 
were presented at the International Fluids conference in Nobeyama, Japan. A report will be 
written and submitted to the conference proceedings. 

Work is still in progress with S. Biringen and F. Hatay (Colorado State University) on the 
direct numerical simulation of compressible Couette flow. During the transition workshop, 
Hatay and I checked the linear stability theory, and he has gone on to do the nonlinear 
calculations. 


James F. Geer 

Work is continuing on hybrid perturbation/variational techniques for solving systems of 
either ordinary or partial differential equations. For ODE’s, we are currently applying the 
methods to several classes of nonlinear vibration problems, with special emphasis on resonant 
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frequency calculations for systems of equations written in Lagrangian or Hamiltonian form. 
For PDF’s, we are investigating how the method might be applied to some exterior boundary 
value problems for elliptic PDF’s, when the boundary of the domain D has an irregular 
shape. The basic idea is to treat D as a perturbation of a simple domain R (e.g. a circle 
or an ellipse, in two dimensions) by embedding it in a one parameter family of domains 
indexed by a parameter a, where 0 < a < 1. Currently, we are applying the method to some 
eigenvalue problems for irregularly shaped domains and to some problems of flows about 
geometrically complicated bodies. In addition, we are combining our method with some 
homotopy methods to further enhance the accuracy of the results we obtain. In particular, 
we hope this will lead to a new extension of slender body theory. So far, the preliminary 
results have been encouraging. Work on the method itself is being done with C. Andersen 
(College of William and Mary), while possible applications are being discussed with E. Liu 
(Fluid Mechanics Division, LaRC). 

The problem of describing aerodynamically generated sound from compact sources of 
vorticity is being studied from a perturbation point of view, with the eddy Mach number 
appearing as the perturbation parameter. The problem appears to be well suited to an ap- 
plication of a slightly modified version of the method of multiple scales. In particular, the 
problem of determining the sound generated by a sphere executing arbitrary (small) oscilla- 
tions in zero mean flow conditions has been investigated successfully by this technique. Once 
the perturbation solution is fully understood, the hybrid perturbation Galerkin technique 
will be applied to it, with the goal of extending the usefulness (and accuracy) of the solution 
to higher Mach numbers than is possible using the perturbation solution alone. This work 
is being carried out with J. Hardin (Acoustics Division, LaRC). 


David Gottlieb 

We have continued our work on the asymptotic stability of high order finite difference 
schemes. In an ICASE report (together with Abarbanel and Carpenter )we reported the new 
phenomenon that high order schemes tend to allow solutions that grow in time. The reason 
is that the inflow boundary conditions are more involved. We try to solve this problem in 
several ways. The first approach is to use penalty type boundary conditions to control this 
growth in time. Using this approach we are testing now a compact fourth order scheme 
with penalty boundary conditions that seems to be asymptotically stable even for systems of 
equations. The second approach is to construct boundary conditions in such a way that the 
scheme satisfies the same energy estimate as the differential equation itself. In joint work 
with Gustafsson we showed that it can be done for explicit schemes but not for compact 
schemes. We feel that the problem has to be answered before compact schemes can be used 
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for problems involving long time integration. 

We have solved the classical problem of the Gibbs phenomenon. It is well known that if 
we are given the first N Fourier modes of an analytic and periodic function then the Fourier 
approximation, based on modes converges exponentially. Gibbs (in 1898) showed that if 
the function is not periodic, then the rate of convergence is ^ away from the discontinuity 
and there is a constant overshoot in the neighborhood of the discontinuity (in this case the 
boundary). 

We have shown that one could construct, with exponential accuracy, the coefficients of 
the expansion of the function based on the Gegenbauer polynomials c^. Thus we can recover 
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the exponential accuracy up to and including the discontinuities. 

Ami Harten 

In collaboration with S. R. Chakravarthy (Rockwell International Science Center) we have 
developed a class of ENO schemes for the numerical solution of multidimensional hyperbolic 
systems of conservation laws in structured and unstructured grids. This is a class of shock- 
capturing schemes which is designed to compute cell-averages to high-order accuracy. The 
ENO scheme is composed of a piecewise-polynomial reconstruction of the solution from 
its given cell- averages, approximate evolution of the resulting initial-value problem, and 
averaging of this approximate solution over each cell. The reconstruction algorithm is based 
on an adaptive selection of the stencil for each cell so as to avoid spurious oscillations near 
discontinuities while achieving high order accuracy away from them. 

I have begun to work on multi-resolution analysis for ENO schemes. Given a function 
u{x) which is represented by its cell-averages in cells which are formed by some unstructured 
grid, we decompose the function into various scales of variation. This is done by considering 
a set of nested grids in which the given grid is the finest, and identifying in each locality the 
coarsest grid in the set from which u(a:) can be recovered to a prescribed accuracy. 

Multi-resolution analysis is applied to ENO schemes in order to reduce the number of 
numerical flux computations which is needed in order to advance the solution by one time- 
step. This is accomplished by decomposing the numerical solution at the beginning of 
each time-step into levels of resolution, and performing the computation in each locality 
at the appropriate coarser grid. We are studying an efficient algorithm for implementing 
this concept in the one- dimensional case; this algorithm can and will be extended to the 
multi-dimensional case with cartesian grids. 
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Thomas L. Jackson 

Work focuses on the time evolution of the interaction of a diffusion flame with a vortex in 
a weak acoustic field, directly applicable to the large scale structures found in compressible 
mixing layers. In particular, the ignition process is under current investigation. In addition, 
the effects of compressibility will be studied. A combination of asymptotics and numerics will 
be used to reduce the complex problem to a model problem, thus isolating key physical effects 
for analysis. This work is in collaboration with M. Macaraeg (Fluid Mechanics Division, 
LaRC) and M.Y. Hussaini. 


Peter A. Jacobs 

A single-block Navier-Stokes solver was completed and was documented in Interim Report 
18 and Report 91-75. The code was then applied to the simulation of transient flow in a 
reflected-shock tunnel with a large Mach number nozzle and results were reported in Report 
91-60. An interesting feature found in the simulation was the formation of a vortex near the 
tube centreline just upstream of the nozzle throat. Theoretical models presently used for the 
prediction of test-gas contamination in reflected-shock tunnels may need to be modified to 
include the effects of the vortex. However, further work is required to determine sensitivity 
of the simulations to grid distortions and facility operating conditions. 

The code was also used to simulate the operation of the Langley expansion tube with 
helium as the working gas. Preliminary results were reported in Interim Report 20. The 
agreement between experimental measurements and the simulation results indicates that we 
now have a reliable tool which may be used to study some of the multidimensional flow 
processes in expansion tubes. An understanding of these processes is key to improving the 
quality of the test flow produced. Since each simulation required some 50 hours of CPU time 
on a Cray-Y /MP and most of the parameter space has yet to be explored, we may attempt 
a parallel implementation of the code in the near future. 


Ashwani Kapila 

Under appropriate conditions, a compressive pulse propagating through an exothermi- 
cally reacting atmosphere steepens and amplifies, eventually resulting in a local thermal 
runaway and a shock. In collaboration with T. Jackson, I have been studying the post- 
runaway evolution of the pulse. Of special interest is the case when the runaway occurs on 
the shock itself. Our analysis is asymptotic, and is a portion of a continuing project on the 
transient events that occur prior to the establishment of detonation waves. 
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David E. Keyes 

The divide-and-conquer paradigm of iterative domain decomposition, or substructuring, 
has become a practical tool in computational fluid dynamics applications because of its flex- 
ibility in accommodating adaptive refinement through locally quasi-uniform grids, its ability 
to exploit multiple discretizations of the operator equations, and the modular pathway it 
provides towards parallelism. Two-scale domain-decomposed preconditionings consisting of 
a coarse grid of subdomain vertices and a large number of independent subdomain problems 
are strong enough to provide convergence rates within a log-factor of optimal as the grid is 
refined, without requiring data exchanges between subdomains on a full hierarchy of scales at 
each iteration. Thus, such algorithms are promising for loosely coupled distributed memory 
computers with high communication latency, including workstation networks. 

Recently we have been seeking generalizations of the diffusively-dominated work de- 
scribed in earlier ICASE Reports 91-19 and 91-20 to high Reynolds number conditions and 
indefinite (Helmholtz) problems. 


Marten T. Landahl 

Work has primarily been directed at analysing the role of algebraic instability on streak 
formation and flow structure in the near-wall region of a turbulent boundary layer. To this 
purpose I have looked at the space-time evolution of a localized three-dimensional disturbance 
in a strong shear flow, paying particular attention to viscous and nonlinear effects. 


D. Glenn Lasseigne 

New research involves the response of a boundary layer behind a nonreacting shock wave 
supported by a wedge. Disturbances are produced as the wedge shape changes or as vorticity 
disturbances in the mean field are propagated through the shock. 

Continuing research focuses on the interactions of disturbances in the flow field with both 
reacting and non-reacting shocks. A combination of asymptotics and numerics is used to 
reduce complex problems to model problems, thus isolating key physical effects for analysis. 
Specific problems include: the non-linear interactions of vorticity /detonation waves, the 
effects of heat release on the response of an oblique detonation subjected to wedge oscillations, 
and the coupling of upstream disturbances and wedge oscillations on the stability of an 
oblique detonation. The goal is to develop numerical routines which capture the instability 
modes predicted by linear theory. Other problems include the effects of streamwise vorticity 
on the stagnation point flow and the effects of more general disturbances on a flame in the 
viscous stagnation region. 
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Dimitri Mavriplis 

Work is continuing on the use of unstructured meshes for solving computational fluid 
dynamics problems in both two and three dimensions. 

In two dimensions, an implicit solver for the full Navier Stokes equations on unstructured 
meshes has been developed in conjunction with V. Venkatakrishnan (Computer Science Cor- 
poration). The method consists of an ILU (incomplete LU factorization) preconditioned 
GMRES implicit iterative technique. In a recent paper (ICASE Report 91-40) this strategy 
was found to provide superior efficiency compared to other preconditioners and other implicit 
iteration techniques, while approaching the computational efficiency of previously developed 
unstructured multigrid techniques. 

In three dimensions, an inviscid flow solver has been developed incorporating adaptive 
meshing techniques and an unstructured multigrid algorithm, thus enabling the efficient and 
accurate computation of three-dimensional flow fields about complex configurations (ICASE 
Report 91-40). Work is continuing on the demonstration of the efficiency of this technique 
with the application of a large scale computation about a complete aircraft configuration 
using a fine mesh of over 0.6 million points. 

Future work will be directed at evaluating new turbulence models for the two-dimensional 
problems. In three dimensions, the investigation of new more efficient and robust gridding 
and adaptive mesh refinement techniques will be pursued. 

In conjunction with J. Saltz, Subhendu Das, and Ravi Ponnusamy, the three-dimensional 
unstructured multigrid Euler solver (non-adaptive) has been implemented on the Intel Hy- 
percube using the PARTI primitives. This work is aimed at demonstrating the applicability 
of the PARTI primitives to practical implementations, and also to employ experience gained 
from particular implementations to aid in the design or modification of the primitives. A 
large scale computation over a complete aircraft configuration running on the iPSC 860 
is in the planning stages in order to demonstrate the efficiency of the parallelized code. 
The implementation will be extended in the future to include the use of adaptive meshing 
strategies. 


Piyush Mehrotra 

Current research in languages and compilers for distributed memory machines has fo- 
cussed on exclusively exploiting data parallelism present in algorithms. However, in addition 
to this loop level parallelism, scientific applications also exhibit functional parallelism at a 
coarser level. We have been investigating an integrated approach which allows the user to 
efficiently exploit both levels of parallelism on distributed memory architectures. In pursuing 
this approach, several issues need to be resolved both at the user (source) level and also at 
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the system/runtime level. 

At the language level, we need to define models/constructs which allow the user to specify 
both the coarse grain and the fine grain parallelism so that the two do not conflict with each 
other. The main objective here is to provide minimal extensions to FORTRAN which strike 
a balance between allowing the user to control factors which critically affect performance 
and hiding the low level details which can be automatically handled by the compiler/runtime 
system. At the system level the most critical issue is the management of resources such as 
processors and memory. In the above model, we have loosely synchronous groups of processes 
interacting asynchronously with other such groups. Each group of processes needs to be 
scheduled across a set of processors simultaneously both from the point of avoiding deadlock 
and also from the point of efficiency. Messages between such groups involve distributed data; 
thus massaging from one distribution to another before the messages are delivered may be 
required. 

The overall objective of the project is to provide a system which allows all levels of paral- 
lelism in a scientific application to be efficiently mapped onto distributed memory machines. 

In joint effort with K.-Y. Wang (Purdue University), we have been investigating an 
optimization technique called message consolidation which attempts to reduce the commu- 
nication and synchronization cost by decreasing the number of messages sent through the 
communication network. However, in certain situations, such combining of messages may 
increase the data synchronization delays thus leading to poorer overall performance. We 
have designed an algorithm for message consolidation and have been investigating heuristics 
for recognizing situations where such a transformation can be profitably applied. 


H. S. Mukunda 

An opposed jet diffusion flame (OJDF) model with finite rate chemistry is being used to 
explain the weak thermal development of a reacting high speed hydrogen-air mixing layer 
which shows small departure of the reaction rates from equilibrium. The use of modified 
flame surface theory permits the reduction to the solution of a set of algebraic equations 
involving chemistry of seven species and nine reversible reactions. The use of this simple 
theory is justified by noting that numerical calculations with full chemistry show a single 
sharp flame front. The solution of the equations is first compared with recent data on 
hydrogen-air OJDF obtained at LaRC. The stretch parameters are extracted from the direct 
numerical simulations of the reacting mixing layer with full chemistry. The results of OJDF 
for these parameters clearly indicate that even though peak temperatures and species mass 
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fractions depart from adiabatic values very significantly, the departures of the reaction rates 
from equilibrium are not large. The effects of Lewis numbers on the flame behavior at high 
temperatures are extracted from the calculations. 


Naomi Naik and John Van Rosendale 

We are concentrating on understanding boundary layer issues in Navier Stokes codes. 
We analyzed and experimentally tested various multiple semi-coarse grid methods, both 
sequential and parallel versions, for their suitability for thin cells and stretched grids. We are 
also completing some past work on spectral preconditioning and determining the complexity 
of concurrent multigrid which will appear as ICASE reports. 

Work on robust parallel multigrid algorithms is continuing in collaboration with J. Dendy 
(Los Alamos National Laboratory). As is well known, to achieve fast convergence on general 
tensor product grids, multigrid algorithms based on line or plane relaxation must be used. In 
three dimensions, only plane-relaxation suffices, in general. Unfortunately, plane relaxation 
multigrid is very difficult to parallelize effectively. Several groups have designed experimental 
parallel implementations of such algorithms, with mixed results. 

A better alternative for achieving robust parallel elliptic solvers is to use multiple coarse 
grids, formed by semi-coarsening the grid separately in each of the coordinate directions, 
together with a simple point-relaxation smoother. This idea, closely related to the hyperbolic 
algorithm of W. Mulder, can be used to create robust and effective elliptic solvers. A class 
of such algorithms now exists, differing mainly in the ways the separate coarse grid solutions 


are combined. 


In recent work, we have established the rapid convergence of the first of these algorithms, 
in which the solutions from separate coarse grids are combined via a local “switch,” based 


on the strength of the discrete operator in each coordinate direction. An improved version 
of this algorithm now exists, free of the necessity of computing the local switch, but we have 


not yet established its rapid convergence. 

In related work, we are developing a variant of this algorithm for problems having severe 
coefficient jumps and strong grid stretching. By combining the new algorithm with inter- 
polation formulas developed previously for plane relaxation algorithms, we can now handle 
problems in which the coefficients vary by a factor of a thousand, from cell to cell, and can 
also handle exponentially stretched grids. 

A variant of this code is being implemented on the Thinking Machines CM 2 at Los 
Alamos National Laboratory. While this class of multigrid algorithms is relatively communi- 
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cation intensive, it seems to map very nicely to SIMD parallel architectures. A hypercube 
version of this algorithm is also being considered, though the mapping issues involved are 
sufficiently different that plane relaxation algorithms may be superior on a hypercube. 


David Nicol 

We are studying “template” communication routing policies for SIMD parallelized nu- 
merical calculations with irregular communication patterns. We have identified conditions 
on the templates under which any work-conserving routing scheme is optimal, and have 
developed an algorithm which is optimal for a constrained class of templates. 

In collaboration with A. Greenberg and B. Lubachevsky (AT&T Bell Labs), algorithms 
for executing trace-driven cache simulations on SIMD architectures are under study. We 
have developed algorithms suitable for simulating a general class of stack replacement policy 
algorithms, as well as random replacement. Previous results have been extended to include 
the efficient analysis of caches with very large set sizes. 

In another effort, we are developing parallel MIMD and SIMD algorithms for the discrete- 
event simulation of networks where messages may be re-routed in the presence of network 
congestion. 

Together with R. Simha (College of William and Mary) and D. Towsley (University of 
Massachusetts), we are looking at the application of the theory of majorization to problems 
in statically load balancing stochastically evolving workloads. Majorization permits us to 
assert that one mapping is better than another in a variety of ways, including expected 
execution time, variation in execution time, space-time product, and reliability. 


Stephen Otto 

The main concern of this study is a greater understanding of vortex wave interactions, 
particularly for situations involving the flow over a plate of variable curvature. We con- 
sider the imposition of a pressure gradient and the effect it has on the onset, and nonlinear 
development of a Gortler vortex mechanism. In joint work with P. Hall, we consider both 
incompressible and compressible boundary layers. We are interested primarily in finite am- 
plitude vortices with small wavelengths, similar to those discussed in Hall & Lakin (1988). 
Also of considerable interest is the possibility that these flow situations may be susceptible to 
mviscid instabilities. For neutral Rayleigh waves to exist we require the flow field to include 
an inflection point. In the incompressible case this is only true if the pressure gradient takes 
a certain form. For a compressible flow situation we find that for high enough Mach numbers 
the flow may contain a generalised inflection point. It is possible to extend this analysis so 
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that the Rayleigh wave becomes large enough to alter the position of the vortex, although 
the wave itself remains linear in character. There is also interest in the verification of these 
results using fully numerical methods. 

A weakly nonlinear study is being made on the effect of crossflow on the most dangerous 
Gortler mode, with a desire to extend this analysis to the fully nonlinear regime; this work 
is joint with A. P. Bassom (University of Exeter, UK). 


Yuh-Roung Ou 

In collaboration with J. Burns, the work concerning active control of unsteady flow by 
moving boundary mechanisms is continuing. An algorithm for computing the viscous flow 
past a rotating cylinder was developed for the purpose of studying the control of cylinder 
forces. Several fundamental types of rotation are demonstrated to be intimately relevant 
to the forces on the cylinder surface. Under a particular type of time-periodic rotation, a 
considerable improvement of the forces has been demonstrated. The computational results 
verify the significance of the proper choice of the rotation rate. Very precise periodicity of 
the force evolutions at certain cases are established, and it may provide direct implication 
for controlling the vortex formation in the cylinder wake. For the case of constant speed of 
rotation, two optimal control problems were considered and solved computationally. 

By treating the rotation rate as a control variable in this model, we will eventually be 
interested in finding the optimal control (i.e. the optimal trajectory of the rotation rate) that 
mi.ximizes the lift-to-drag ratio over a fixed time interval. Although in this study, an optimal 
control problem associated with the constant rotation rate was solved by direct computation, 
it is still important to explore further the implementation of a computational algorithm to 
calculate the optimal solution from the full necessary conditions. This approach may be 
applied to practical control problems arising in several areas of fluid flow. 


Terry Pratt 

A promising new method, called the KCD method, for automatic generation of an MIMD, 
distributed memory, parallel program from existing sequential code is being studied. The 
method is based on a theoretical model called ’’kernel-control decomposition” (KCD). The 
KCD method uses a source-to-source translation of the sequential code into a ’’control” 
process and multiple ’’kernel” processes. The control process determines the control path 
and assigns work (in the form of ’’kernel segments”) to the kernel processes, which operate 

in parallel. 

The KCD method has been applied to three applications codes to date, the NASA 
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LAURA code (Langley Aerothermodynamic Upwind Relaxation Algorithm), developed by 
R GnofFo (Space Systems Division, LaRC) and two of the codes distributed as part of the 
Perfect Club benchmark suite (FLO 52 and TRFD). Measurements of the performance of the 
parallel versions of these codes are being made on the 32-node NASA/ICASE Intel iPSC/860. 


Peter Protzel 

Together with D. Palumbo (Information Systems Division, LaRC) we have completed a 
study to investigate the potential benefits of using Artificial Neural Networks (ANNs) as 
components of critical, real time systems with high reliability requirements. The first part 
of this study considered an example of an application in which an ANN controlled the re- 
allocation of tasks in a fault-tolerant, distributed multiprocessor system. This reallocation 
is necessary if one of the processors fails and all tasks have to be distributed among the 
remaining processors as quickly as possible by observing certain constraints and by approx- 
imately balancing the load of the processors. Since the ANN becomes a critical component 
of the system in this application, we studied the fault-tolerance of the ANN by simulated 
fault-injection experiments. The results show a surprising degree of fault- tolerance with a 
graceful performance degradation even after multiple faults. Thus, an ANN implemented 
as an analog VLSI chip is an attractive alternative for performing certain critical functions, 
especially on long-term missions where cumulative component failures have to be expected 
without the possibility of repair or maintenance. More details on this application can be 
found in ICASE Report No. 91-45 and in a NASA Technical Paper currently in print. 

The second part of our project focused on a different type of ANN to control the posi- 
tioning of a robot arm with five degrees of freedom and with visual feedback provided by 
two cameras. This example application was chosen to demonstrate the fault-tolerance and 
adaptivity that can be achieved by using an ANN as the controller. The network learns and 
adapts on-line without any outside intervention. The continuous feedback provided by the 
cameras is used to generate an average positioning error that the networks tries to minimize. 
Three different issues were addressed from a fault-tolerant systems perspective; a possible 
(partial) failure of the controller, a failure of the robot arm, and a sudden change in the 
environmental or external conditions. In our simulations we showed, for example, how the 
system recovers from sudden failures of one or two arm joints, from a loss of 30% of the 
network’s neurons, and from geometric changes in the camera positions. These remarkable 
results lead us to conclude that autonomous systems represent one of the most promising 
application areas for this kind of self-supervising ANN. Except for the fault-tolerance of the 
ANN itself, the main benefits can be achieved even if the ANN is only simulated and treated 
as an algorithm which can be executed by a standard microprocessor. 


18 


Janies Quirk 

Work continued on the development of an adaptive grid algorithm for simulating shock 
hydrodynamic phenomena. Previously it was shown that the algorithm required only mod- 
est computing resources in order to produce very high resolution simulations. The basic 
efficiency of the algorithm stems from the fact that it refines in time as well as space. A 
larger number of smaller time steps are taken for fine mesh cells than for coarse mesh cells. 
Thus the scheme is inherently well suited to the task of producing time accurate simulations 
for flows which contain disparate length scales. However, a certain amount of geometric 
flexibility had been sacrificed so that the spatial and temporal refinement strategies could 
be linked harmoniously. Recent work has concentrated on removing this limitation. 

Whereas in the past the algorithm essentially employed a body-fitted grid, it now uses 
a grid which intersects solid surfaces. In essence, a solid body is allowed to blank out areas 
of a cartesian grid. This results in a number of cut cells along the boundary of the solid 
which are meted out for special treatment during the integration of the flow solution. In 
principle such a strategy can cope with arbitrarily complicated bodies. Initial results are 
encouraging; at least for simple shapes such as wedges and cylinders the cut-cell approach 
produces results that are comparable to those produced by the body-fitted code. Further 
work is required to test this new approach on the sorts of complicated shapes that can be 
tackled routinely using unstructured grids. 

Looking to the future, we intend to investigate the extent to which our adaptive grid 
algorithm can exploit the power offered by parallel computers, for this extent will have an 
important bearing on the long term usefulness of the scheme. 


Rolf Radespiel and John Van Rosendale 

Most flow solvers used in current aerospace programs exhibit slow convergence to the 
steady-state solution, resulting in high computer costs, large turn-around times, and slow 
progress towards the desired aerospace vehicle design. The reason is the appearance of 
flow phenomena with very different scales, also complicated by nonlinear behavior. Over 
the past several years, Radespiel has been developing a class of multigrid algorithms, using 
multiple semicoarsened grids, for the hypersonic regime. While the multiple coarse grids 
required are more expensive, such multigrid solvers appear much faster in the context of the 
Navier Stokes equations than traditional (full coarsening) multigrid solvers. The reason for 
this improvement seems to be that the multiple coarse grids effectively remove the stiffness 
arising from high aspect ratio cells in the boundary layer. 

In order to understand such phenomena, we have been applying Fourier analysis to a 
two-level multigrid algorithm for the two-dimensional convection equation in collaboration 
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with C. Swanson (Fluid Mechanics Division, LaRC). This is yielding insight, though there j 

are a number of complicating effects, including flow alignment issues, and the problem of | 

the different speeds of convective and acoustic waves. The Runge Kutta “smoothers,” used 
in Navier Stokes codes, are also difficult to analyze, since the smoother and the multigrid 
algorithm interact in complicated ways. 

The two dimensional semicoarsening multigrid algorithms for the Navier Stokes equations 
developed so far are not directly usable in three dimensions, for several reasons. The primary 
reason is the great cost, in both storage and computation time, of all of the coarse grids 
required in three dimensions. In the elliptic case, a simple solution is to coarsen by a factor 
of four rather than two in each coordinate direction. Coarsening by a factor of two is 
expensive, since the multigrid algorithm requires eight times the amount of storage of the 
fine grid, and proportionately more work as well. Coarsening by a factor of four implies 
storage only about double that of the fine grid, but requires a better smoothing iteration, 
since more Fourier modes need to be effectively damped by the smoother. 

In coming months we plan to implement a multigrid algorithm using the multiple coarse 
grids produced by factor of four coarsening in each coordinate direction, first for the elliptic 
case, then for the full Navier Stokes equations. In the elliptic case, we now have several 
effective ways of combining the corrections from the different coarse grids, all performing 
about equally well. This is not at all true in the hyperbolic case, where the way in which the 
results from the separate coarse grids are combined appears to be critical. The numerical 
test codes currently in progress should resolve the issue of how best to do this. 


Dan Reed and David Nicol 

We are developing simulation models and analytic models of dynamic load-balancing 
policies that choose a destination processor at random. Such policies are easy to implement, 
although they ignore state information. Our aim is to assess the quality and costs of these 
types of policies. 


Philip Roe 

I concentrated on clarifying the scalar update algorithms for use in the new class of 
multidimensional upwind schemes. This led to a unification of previous work and to some 
estimates of truncation error for both the linear and nonlinear cases. 

I also did some preliminary work to explore the use of CFD methods for aeroacoustic 
problems, and looked at a promising technique for solving the steady Navier-Stokes equations 
by fast explicit methods. 
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Joel Saltz, Scott Berryman, Raja Das, Dimitri Mavriplis, Vijay Naik, Ravi Pon- 
nusamy and Sesh Venugopal 

We have further developed and refined Fortran and C callable distributed memory PARTI 
procedures to carry out execution time preprocessing needed for sparse matrix and unstruc- 
tured mesh problems. Our new optimizations include methods for increasing the number of 
situations in which it is simple to avoid communication by reusing stored off-processor data. 
We also have developed a suite of primitives that track runtime data dependencies. These 
primitives are useful in such problems as sparse direct methods. We have used our new 
primitives to implement on the iPSC/860 applications such as Mavriplis’ 3-D unstructured 
mesh explicit and multigrid Euler solver, Horst Simon’s recursive spectral domain parti- 
tioner, a conjugate gradient solver, a substantial portion of the molecular dynamics program 
CHARMM and a sparse direct systems solver. 

The explicit Euler solver on a 210K mesh carries out calculations at a rate in excess 
of 160 Mfiops on 64 processors. The multigrid Euler solver carries out calculations at a 
rate in excess of 130 Mfiops on 64 processors. Most of the time in both calculations is 
spent on computation. We have promising results that indicate that reordering the data and 
loop iterations will lead to a significant improvement in computational rates. The explicit 
unstructured Euler solver has been ported to the Delta machine at CalTech; performance 
studies are under way. 

The primitives are available now from netlib or by anonymous ftp from ra.cs.yale.edu. A 
manual for some primitives is available as ICASE Interim report 91-17, other manuals are 
pending. 


Joel Saltz, Craig Chase and Kay Crowley 

Fortran and C callable procedures which carry out symbolic and execution time prepro- 
cessing needed for data exchange in block structured and structured adaptive programs are 
under study. We have developed a new suite of primitives for block structured and structured 
adaptive problems and have tested the primitives using routines from the block structured 
Navier Stokes code written by V. Vatsa (Fluid Mechanics Division, LaRC). We have also de- 
veloped a rudimentary compiler capable of generating code that runs on distributed memory 
processors. This compiler has been tested for simple block structured kernels. 

Earlier versions of the primitives were used in a structured adaptive PDE solver developed 
by J. Scroggs; currently the primitives are being used to implement a NASA block structured 
Euler solver. The primitives are being used to implement and extend Fortran D in joint work 
with the groups of Ken Kennedy and Geoffrey Fox at the NSF Center for Research on Parallel 
Computing. 
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Joel Saltz, Raja Das and Ravi Ponnusamy 

We have developed two methods which we believe will play an important role in any 
distributed memory compiler able to handle sparse and unstructured problems. We described 
in ICASE report 91-73 how to link runtime partitioners to distributed memory compilers. 
In our scheme, programmers can implicitly specify how data and loop iterations are to be 
distributed between processors. This insulates users from having to deal explicitly with 
potentially complex algorithms that carry out work and data partitioning. 

In ICASE report 91-73 we also described a viable mechanism for tracking and reusing 
copies of off-processor data. In many programs, several loops access the same off-processor 
memory locations. As long as it can be verified that the values assigned to off- processor 
memory locations remain unmodified, we showed that we can effectively reuse stored off- 
processor data. We presented experimental data from a 3-D unstructured Euler solver run 
on an iPSC/860 to demonstrate the usefulness of our methods. 

The Parti primitives can be used directly by programmers or incorporated into compiler 
for distributed memory architectures. They are currently being used to implement and 
extend Fortran D in joint work with K. Kennedy, C. Koelbel and R. von Hanxleden (Rice 
University) and G. Fox and A. Choudhary (Syracuse University). 


Joel Saltz and Satyanarayan Gupta 

We have developed PARTI primitives for the Maspar and have used these primitives to 
implement Dimitri Mavriplis’ 3-D explicit unstructured Euler code. The Parti code must 
preprocess message patterns in a manner that allows the Maspar router to be effectively used. 
We are currently investigating heuristics that can be used to carry out this preprocessing. 


Sutanu Sarkar 

We are engaged in the direct simulation and Reynolds stress modeling of turbulent flows. 
Our simulations of homogeneous shear flow show that compressibility leads to a significant 
decrease in the growth rate of turbulent kinetic energy. Augmented dissipation (through the 
dilatational dissipation) and energy transfer to pressure fluctuations (through the pressure- 
dilatation) have been identified as mechanisms responsible for the decreased growth rate. 
Furthermore, we have proposed an algebraic model for the pressure-dilatation to be used in 
Reynolds stress modeling of high-speed flows. To complement the non-linear direct simula- 
tions, we have also studied the linear stability of compressible homogeneous shear flow. 
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Compressibility effects on turbulence structure are now being obtained in collaboration with 
G. Erlebacher by analysis of the shear flow simulations. We are also addressing issues related 
to turbulence anisotropy and pressure-strain modeling for compressible shear flows. 


Jeffrey S. Scroggs and Peter A. Jacobs 

An algorithm is being developed for the driver of a a multi-domain PDE solver, with 
target application being the time-accurate simulation of high-speed gasdynamics using struc- 
tured grids for each of the domains. The driver algorithm uses a time-accurate solver for 
a single block essentially as a black-box. In addition to the physical boundary conditions, 
multidomain solvers require sharing information at the (artificial) interfaces between do- 
mains. We allow different time steps to be taken on adjacent grids. When Atij Atj > 1, we 
interpolate the solution from grid j to grid i, and will use a flux-correction scheme to share 
data from grid i to grid j . The interpolation has been implemented, and the implementation 
of the flux-correction is a subject of current research. 

We demonstrated multi-driver for a three-block decomposition of axi-symmetric flow 
around the leading edge of a cone using the single-block integrator described in previous 
reports by P. Jacobs. 


Sharon O. Seddougui 

The instability of hypersonic flow pcist a cone is being investigated with A. Bassom 
(University of Exeter). The effect of an attached shock is studied in the weak interaction 
region. We find that the triple deck structure is identical to that given by Cowley and 
Hall (1990) for hypersonic flow over a wedge. In this case the local position on the cone is 
determined so that the effects of curvature are important. Solutions of the dispersion relation 
are currently being obtained. Joint work with A. Bassom was completed investigating the 
effect of suction on the stability of the viscous stationary modes for the flow over a rotating 
disc. MacKerrell (1988) showed that this crossflow instability is subcriticaJ, and that there 
exists a threshold value of the disturbance amplitude. We show that suction does not alter 
these results although it lowers the value of the threshold amplitude. Thus suction increases 
the likelihood of this mode in an experimental setting. Asymptotic results for large suction 
and large blowing were obtained. The extension of this problem into the compressible regime 
is currently being studied. 
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Charles G. Speziale 

An analysis of the Renormalization Group (RNG) K — e model of Yakhot & Orszag 
was conducted in collaboration with T.B. Gatski (Fluid Mechanics Division, LaRC) and S. 
Thangam. It was found that with a correction of the numerical value of one of the constants 
- and with the incorporation of a nonlinear strain dependent modification to the production 
of dissipation term - a substantially better description of turbulent shear flows was achieved. 
In particular, excellent predictions were obtained for homogeneous shear flow and turbulent 
flow over a backward facing step: two test cases that are commonly used to benchmark 
the performance of turbulence models. Future research Is planned on the development of 
further modifications that will allow this two-equation model to fully accommodate the rapid 
distortion limit. 

Research was conducted in collaboration with P.S. Bernard (University of Maryland) on 
the theory of self-preservation in homogeneous turbulence. It was found that self-preservation 
renders an alternative equilibrium structure which includes a asymptotic power law decay 
in isotropic turbulence and a production-equals-dissipation equilibrium in homogeneous shear 
flow. The ability of the self preserving model to accommodate the limit of zero viscosity - 
wherein it predicts a finite time enstrophy blow up - allows for a better description of non- 
equilibrium turbulent flows. Applications of these results to the development of improved 
turbulence models are currently under investigation. 

A new near- wall two-equation turbulence model for supersonic flows has been developed 
in collaboration with R.M.C. So (Arizona State University). A variable density extension of 
the asymptotically consistent near- wall model of So and co-workers for incompressible flows 
is combined with some new dilatational models developed recently by S. Sarkar. Excellent 
results are obtained in the supersonic flat plate turbulent boundary layer for Mach numbers 
as large as 10. Future research is planned where this approach will be extended to second- 
order closure models. 


Shlomo Ta’asan 

The development of efficient multigrid solvers for constraint optimization problems gov- 
erned by partial differential equations has continued with research in two directions. The 
first, which is well developed by now, deals with problems in which the parameter space 
on which optimization is done is of finite dimension in the differential formulation of the 
problem. The methods use relaxation for the parameter space in a multilevel way. Param- 
eters that have a non-smooth effect on the solution are relaxed on fine levels while those 
of smooth effect are solved for on coarse grids only. The methods use adjoint variables to 
define a descent direction for the minimization problem. The other direction focuses on 
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problems in which the optimization is over an infinite dimensional parameter space (in the 
differential level). Also here the same type of ideas for the treatment of the different scales 
in the problems is being used. Experiments with some model problems involving elliptic 
partial differential equations as the constraint equations have been performed showing that 
the full optimization problem can be solved with a computational cost which is only a few 
times more than that of solving the PDE alone. 

The above ideas have been applied in aerodynamics design problems where airfoils are 
to be calculated so as to meet certain design requirements, for example, to give pressure 
distribution in some flow conditions which are closest to a given pressure distribution. The 
present model for the flow is the transonic small disturbance equation in which an airfoil 
is modelled by the small disturbance boundary condition. The shape of the airfoil in these 
calculations is being expressed in terms of a finite number of given shape functions with 
amplitudes to be found by the design process. We have shown that subsonic design problems 
can be solved in a computational cost which is just a few times (2-3) that of the flow solver. 
This work is jointly done with M. D. Salas (Fluid Dynamics Division, LaRC) and G. Kuruvila 
(Vigyan, Inc.). 

New multigrid solvers for inviscid flow problems are being developed in which the conver- 
gence rates are to be like that of the full potential equation. These convergence rates have 
been predicted by Fourier analysis. These solvers employ relajcation methods of the Gauss- 
Seidel type with a proper modification to handle systems of partial differential equations. 
At present we focus on the incompressible inviscid case, working with body fitted grids in 
twc dimension, using MAC discretization schemes. 

Another area of research is the development of new computational techniques in elastic- 
plastic problems. Here a set of non-linear equations which involve also weak dependence of 
the history of the loading are to be calculated. A new numerical scheme for the integration 
of the plastic forces has been constructed. Using this scheme an efficient multigrid solver 
is under investigation. It uses appropriate fine to coarse grid transfer of the yield functions 
so as to enable the solution of the large scale behavior of the plastic forces to be done on 
the coarsest levels, yielding the accuracy of the fine grid solution. The method is under 
study with some real elastic-plastic problems related to the understanding of fatigue failure 
in metals. The convergence rates that are obtained with the present formulation of the 
algorithm are not satisfactory, i.e., much slower than the linear case. Further investigation 
will focus on understanding this point. 


Siva Thangam 

Turbulent separated flow past a backward-facing step is considered from modeling, com- 
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putational and experimental point of view. The first and the second phases of this collabora- 
tive effort involving C. Speziale, M.Y. Hussaini and S. Kjelgaard (Fluid Mechanics Division, 
LaRC) have been completed. In the first phase, the equa,tions of motion based on the two- 
equation Reynolds stress closure model are solved using a finite-volume algorithm. The 
results show that when the boundary conditions are consistently applied, the experimental 
and computational results are in very good agreement. During the second phase, several 
refinements were made to the two-equation (k-e) model to improve its predictive capability. 
It was ascertained that the optimal selection of the model constants, inclusion of anisotropy, 
or the accurate representation of the turbulence production-dissipation balance, can con- 
siderably improve the two-equation model predictions. Computations were also performed 
using the large-eddy simulation method with an eddy-viscosity based subgrid scale model. 
The predictions of the time-averaged results of the large-eddy simulation (with 96x48x32 
mesh) were compared with the available experimental findings as well as with those from 
the two-equation model computations (based on a 100x50 mesh) for a step to channel height 
ratio of 1:3. It was found that the dominant features of the flow field - namely, the size of 
the separated flow region and the mean velocity profiles - can be accurately predicted. In 
addition, during this phase the construction and testing of the backward-facing step facil- 
ity at the Experimental Methods Branch was completed. During the next phase, detailed 
experiments as well as additional computations based on full Reynolds-stress closure model 
and large-eddy simulations will be performed for a variety of flow conditions. The issues 
involving the the subgrid scale representation for the large-eddy simulations as well as the 
predictive capability of the full-Reynolds stress closure model will be addressed. 


Eli Turkel 

Work was continued on extensions of central difference and multigrid techniques to hy- 
personic problems. The standard method does not converge due to the presence of strong 
shocks. By revising the switch near shocks and using residual smoothing on the way to finer 
meshes, converged solutions with good convergence rates were obtained for Mach 10 and 
Mach 20 turbulent flows about many aerodynamic bodies including the HSCT project. A 
sequence of ICASE reports were written for various applications. 

A new project was begun in constructing new implicit schemes for spectral methods that 
are easier to invert than a full implicit trapezoidal rule or backward Euler. These methods are 
only first order accurate for time dependent problems but spectrally accurate for steady state 
problems. In some simple test cases on the one dimensional Euler equations the algorithm 
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converged to the steady state within four time steps on a mesh with 64 nodes using the 
pseudospectral Chebyshev method. The algorithm is currently being extended to two space 
dimensions. 


Bram van Leer 

The work on a multi-dimensional approximate Riemann solver with C. Rumsey (Fluid 
Mechanics Division, LaRC) was completed with Rumsey’s doctoral thesis at the University of 
Michigan. Although some very good results and insights have been obtained, it is felt that in 
this work the amount of information extracted from the two" input states is too large to lead 
to a robust method. In particular, the use of the input data to define a rotated coordinate 
frame in which the Riemann problem based on the same input data is then solved, leads to 
noisy, unsteady solutions where a steady solution should be found. Although one still can 
think of some modifications that are worth trying, these will presently not be implemented. 

The work on local preconditioning of the Euler equations continued in collaboration with 
E. Turkel, and led to quite a surprise: one matrix of a family derived earlier by Turkel has 
exactly the same effect on the wave speeds as the preconditioning matrix derived by Van 
Leer, Lee and Roe. The matrices are different: Van Leer’s is symmetric for the symmetrized 
Euler equations, while Turkel’s is not. There probably is a whole two-parameter family of 
matrices that can achieve the optimal characteristic condition number; an explicit formula 
for these has not yet been obtained. Further work on the preconditioning matrix for stretched 
grids is underway. It is possible to overcome, at lecist in part, the stiffness caused by the 
presence of strongly elongated cells, such as found in boundary layers. While the use of 
the preconditioning for the Navier-Stokes equations still has to be explored, the matter of 
the large cell aspect-ratio can be studied independently with Euler schemes used on Navier- 
Stokes type grids. 


Sesh Venugopal 

One of the objectives of the work with J. Saltz and V. Naik, (I.B.M. T.J. Watson Re- 
search Center.) in progress is to develop a new set of efficient primitives for automating the 
communication during the execution of a task dependency directed acyclic graph (DAG). 
These primitives are designed for run-time support of any computation whose communication 
can be expressed as a coarse-grained DAG, with the nodes in the DAG representing linear 
blocks of data. As a test DAG, we have considered the column oriented sparse Cholesky 
factorization algorithm. I have implemented four versions of this algorithm on the iPSC/860 
and have obtained execution timings for each. One of the versions is based on dynamic 
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task scheduling while the other three use a static or pre-determined ordering of tasks on 
processors. One of the static versions uses PARTI primitives for message passing, another 
performs manual communication and the third uses the new DAG primitives. 

A second objective of this work is to determine the effects of partitioning and scheduling of 
sparse factorization on load balancing and communication (ICASE Report No. 91-80). Work 
is continuing towards developing efficient block partitioning and scheduling schemes. We plan 
to implement the block factorization code on the iPSC/860 and compare its performance 
with column based codes, and extend the DAG primitives to handle higher dimensional 
blocks of data. 


Hans Zima, Piyush Mehrotra and Barbara Chapman 

In recent years, a number of distributed memory multiprocessing computers have been 
introduced into the market (e.g. Intel’s iPSC series, NCUBE, and several transputer based 
systems). In contrast to shared memory machines, these architectures are less expensive to 
build and are potentially scalable to a large number of processors. However, the associated 
programming paradigm requires the user to specify complete details of the synchronization 
and the communication of data between processors as dictated by the algorithm. Experience 
has shown that forcing the user to provide such low level details not only makes distributed 
memory programming tedious and error prone but also inhibits experimentation. 

Based on our earlier work on languages and compilers for distributed memory machines, 
we have been investigating machine-independent language extensions for FORTRAN 77, 
called Vienna Fortran (ICASE Report 91-2). These extensions allow the user to write pro- 
grams for distributed memory architectures using global data references only. Since the 
distribution of data to the processors is critical for performance, the extensions described 
here permit the user to explicitly control the mapping of data onto the processors. Vienna 
Fortran supports a wide range of facilities for specifying data distributions, including dis- 
tribution by alignment, where one array is mapped so that it has a fixed relationship with 
another array. It also supports dynamic redistribution of data. Frequently occurring dis- 
tributions can be specified in a simple manner, whereas the full language permits complex 
mappings. The overall aim of the language extensions provided by Vienna Fortran is to 
make the transition from the sequential algorithm to a parallel version as easy as possible, 
without sacrificing performance. 
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REPORTS AND ABSTRACTS 


Banks, H.T., G. Propst, and R.J. Silcox: A comparison of time domain boundary conditions 
for acoustic waves in wave guides. ICASE Report No. 91-27 , August 15, 1991, 25 pages. 
Submitted to Journal of Sounds and Vibration. 


We consider several types of boundary conditions in the context of time domain models for 
acoustic waves. Experiments with four different duct terminations (hardwall, free radiation, 
foam, wedge) were carried out in a wave duct from which reflection coefficients over a wide 
frequency range were measured. These reflection coefficients are used to estimate parameters 
in the time domain boundary conditions and a comparison of the relative merits of the models 
in describing the data is presented. Boundary conditions which yield a good fit of the model 
to the experimental data were found for all duct terminations except the wedge. 


Lasseigne, D. G., T. L. Jackson, and M. Y. Hussaini: Non-linear interaction of a detona- 
tion/vorticity wave. ICASE Report No. 91-34 , April 1, 1991, 27 pages. To appear in Physics 
of Fluids A. 


The interaction of an oblique, overdriven detonation wave with a vorticity disturbance is 
investigated by a direct two-dimensional numerical simulation using a multi-domain, finite- 
difference solution of the compressible Euler equations. The results are compared to those 
of linear theory, which predict that the effect of exothermicity on the interaction is relatively 
small except possibly near a critical angle where linear theory no longer holds. It is found 
that the steady-state computational results whenever obtained in this study agree with 
the results of linear theory. However, for cases with incident angle near the critical angle, 
moderate disturbance amplitudes, and/or sudden transient encounter with a disturbance, 
the effects of exothermicity is more pronounced than predicted by linear theory. Finally, it 
is found that linear theory correctly determines the critical angle. 


Grosch, C. E., and T. L. Jackson: Ignition and structure of a laminar diffusion flame in 
a compressible mixing layer with finite rate chemistry. ICASE Report No. 91-35 , April 8, 
1991, 28 pages. Submitted to Physics of Fluids A. 


In this paper we consider the ignition and structure of a reacting compressible mixing layer 
using finite rate chemistry lying between two streams of reactants with different freestream 
speeds and temperatures. Numerical integration of the governing equations show that the 
structure of the reacting flow can be quite complicated depending on the magnitude of the 
Zeldovich number. In particular, for sufficiently large Zeldovich number, the three regimes 
first described by Linan and Crespo (1976); ie., ignition, deflagration, and diffusion flame, 
occur in supersonic as well as in subsonic flows. An analysis of both the ignition and diffusion 
flame regimes is presented using a combination of large Zeldovich number asymptotics and 
numerics. This allows us to analyze the behavior of these regimes as a function of the 
parameters of the problem. For the ignition regime, a well defined ignition point will always 
exist provided the adiabatic flame temperature is greater than either freestream temperature. 
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One important result is that at supersonic speeds ignition occurs far downstream from the 
plate and, as the flow is accelerated to hypersonic speeds, ignition is exponentially delayed. 
For the diffusion flame regime, the location of the flame changes significantly with changes 
in the equivalence ratio and the Schmidt numbers. 


Thompson, C. P., G. K. Leaf, and J. Van Rosendale: A dynamically adaptive multigrid 
algorithm for the incompressible Navier-Stokes equations - Validation and model problems. 
ICASE Report No. 91-36 , June 27, 1991, 32 pages. Submitted to J. of Applied Numerical 
Mathematics. 


We describe an algorithm for the solution of the laminar, incompressible Navier-Stokes 
equations. The basic algorithm is a multigrid method based on a robust, box-based smooth- 
ing step. Its most important feature is the incorporation of automatic, dynamic mesh re- 
finement. Using an approximation to the local truncation error to control the refinement, 
we use a form of domain decomposition to introduce patches of finer grid wherever they are 
needed to ensure an accurate solution. This refinement strategy is completely local: regions 
that satisfy our tolerance are unmodified, except when they must be refined to maintain 
reasonable mesh ratios. This locality has the important consequence that boundary layers 
and other regions of sharp transition do not (lqsteai(rq mesh points from surrounding regions 
of smooth flow, in contrast to moving mesh strategies where such (lqstealing(rq is inevitable. 

Our algorithm supports generalized simple domains, that is, any domain defined by 
horizontal and vertical lines. This generality is a natural consequence of our domain decom- 
position approach. We base our program on a standard staggered-grid formulation of the 
Navier-Stokes equations for robustness and efficiency. To ensure discrete mass coiiservation, 
we have introduced special grid transfer operators at grid interfaces in the multigrid algo- 
rithm. While these operators complicate the algorithm somewhat, our approach results in 
exact mass conservation and rapid convergence. 

In this paper, we present results for three model problems: the driven-cavity, a backward- 
facing step, and a sudden expansion/contraction. Our approach obtains convergence rates 
that are independent of grid size and that compare favorably with other multigrid algorithms. 
We also compare the accuracy of our results with other benchmark results. 


Speziale, Charles G., Thomas B. Gatski, and Nessan Fitzmaurice: An analysis of RNG based 
turbulence models for homogeneous shear flow. ICASE Report No. 91-37 , April 16, 1991, 12 
pages. To be submitted to Physics of Fluids A, 


In a recent paper [Phys. Fluids A2:1678-1684, 1990], the authors compared the perfor- 
mance of a variety of turbulence models including the K - e model and the second-order 
closure model derived by Yakhot and Orszag based on Renormalization Group (RNG) meth- 
ods. The performance of these RNG models in homogeneous turbulent shear flow was found 
to be quite poor, apparently due to the value of the constant C^i in the modeled dissipation 
rate equation which was substantially lower than its traditional value. However, recently a 
correction has been made in the RNG based calculation of Cei- It is shown herein that with 
the new value of Cei, the performance of the RNG K - e model is substantially improved. 
On the other hand, while the predictions of the revised RNG second-order closure model are 
better, some lingering problems still remain which can be easily remedied by the addition of 
higher-orde r t erm s. 
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Shu, Chi-Wang, Gordon Erlebacher, Thomas A. Zang, David Whittaker, and Sta,nley Os- 
her: High-order ENO schemes applied to two- and three-dimensional compressible flow. 
ICASE Report No. 91-38 , April 30, 1991, 44 pages. Submitted to J. Applied Numerical 
Mathematics. 


High-order essentially non-oscillatory (ENO) finite-difference schemes are applied to the 
two- and three-dimensional compressible Euler and Navier-Stokes equations. Practical issues, 
such as vectorization, efficiency of coding, cost comparison with other numerical methods 
and accuracy degeneracy effects, are discussed. Numerical examples are provided which are 
representative of computational problems of current interest in transition and turbulence 
physics. These require both non-oscillatory shock capturing and high resolution for detailed 
structures in the smooth regions and demonstrate the advantage of ENO schemes. 


Fu, Yibin B., and Philip Hall: Nonlinear development and secondary instability of Gdrtler 
vortices in hypersonic flows. ICASE Report No. 91-39 , May 2, 1991, 47 pages. Submitted 
to Euro. J. of Mechanics. 


In a hypersonic boundary layer over a wall of variable curvature, the region most sus- 
ceptible to Gortler vortices is the temperature adjustment layer over which the basic state 
temperature decreases monotonically to its free stream value (Hall & Fu (1989), Fu, Hall & 
Blackaby (1990)). Except for a special wall curvature distribution, the evolution of Gortler 
vortices trapped in the temperature adjustment layer will in general be strongly affected by 
boundary layer growth through the 0(M^/^) curvature of the basic state, where M is the 
free stream Mach number. Only when the local wavenumber becomes as large as of order 
do nonparallel effects become negligible in the determination of stability properties. 
In the latter case, Gortler vortices will be trapped in a thin layer of 0(c^^^) thickness which 
is embedded in the temperature adjustment layer; here e is the inverse of the local wavenum- 
ber. In this paper, we first present a weakly nonlinear theory in which the initial nonlinear 
development of Gortler vortices in the neighbourhood of the neutral position is investigated 
and two coupled evolution equations are derived. From these two evolution equations we can 
determine whether the vortices are decaying or growing depending on the sign of a constant 
which is related to the wall curvature and the basic state temperature. In the latter case, it 
is found that the mean flow correction becomes as large as the basic state at distances 0(1) 
downstream of the neutral position. Next, we present a fully nonlinear theory concerning 
the further downstream development of these large- amplitude Gortler vortices. It is shown 
that the vortices spread out across the boundary layer. The upper and lower boundaries 
of the region of vortex activity are determined by a free-boundary problem involving the 
boundary layer equations. Finally, the secondary instability of the flow in the transition 
layers located at the upper and lower edges of the the region of vortex activity is considered. 
The superimposed wavy vortex perturbations are spanwise periodic travelling waves which 
are 7t/2 radians out of phase with the fundamental. The dispersion relation is found to be 
determined by solving two coupled differential equations and it is shown that an infinite 
number of neutrally stable modes may exist. 
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Venkatakrishnan, V., and Dimitri J. Mavriplis: Implicit solvers for unstructured meshes. 
ICASE Report No. 91-40, May 6, 1991, 23 pages. 


We develop and test implicit methods for unstructured mesh computations. The ap- 
proximate system which arises from the Newton-linearization evolution operator is solved by 
using the preconditioned GMRES (Generalized Minimum Residual) technique. We investi- 
gate three different preconditioners, namely, the incomplete LU factorization (ILU), block 
diagonal factorization and the symmetric successive over- relaxation (SSOR). The precondi- 
tioners have been optimized to have good vectorization properties. We also study SSOR 
and ILU themselves as iterative schemes. The various methods are compared over a wide 
range of problems. Ordering of the unknowns, which affects the convergence of these sparse 
matrix iterative methods, is also investigated. Results are presented for inviscid and tur- 
bulent viscous calculations on single and multi-element airfoil configurations using globally 
and adaptively generated meshes. 


Mavriplis, D. J.: Three dimensional unstructured multigrid for the Euler equations. ICASE 
Report No. 91-41 , May 7, 1991, 29 pages. Submitted to AIAA Journal. 

The three-dimensional Euler equations are solved on unstructured tetrahedral meshes 
using a multigrid strategy. The driving algorithm consists of an explicit vertex-based finite- 
element scheme, which employs an edge-based data-structure to assemble the residuals. The 
multigrid approach employs a sequence of independently generated coarse and fine meshes to 
accelerate the convergence to steady-state of the fine grid solution. Variables, residuals and 
corrections are passed back and forth between the various grids of the sequence using linear 
interpolation. The addresses and weights for interpolation are determined in a preprocessing 
stage using an efficient graph traversal algorithm. The preprocessing operation is shown to 
require a negligible fraction of the CPU time required by the overall solution procedure, while 
gains in overall solution efficiencies greater than an order of magnitude are demonstrated on 
meshes containing up to 350,000 vertices. Solutions using globally regenerated fine meshes 
as well as adaptively refined meshes are given. 


Sarkar, S.: Modeling the pressure-dilatation correlation. ICASE Report No. 91-42 , May 7, 
1991, 20 pages. Submitted to Physics of Fluids A. 


It is generally accepted that the pressure-dilatation, which is an additional compressibility 
term in the turbulence transport equations, may be important for high-speed flows. Recent 
direct simulations of homogeneous shear turbulence have given concrete evidence that the 
pressure-dilatation is important insofar that it contributes to the reduced growth of turbulent 
kinetic energy due to compressiblity effects. The present work addresses the problem of 
modeling the pressure-dilatation. We first isolate a component of the pressure-dilatation 
which exhibits temporal oscillations and, using direct numerical simulations of homogeneous 
shear turbulence and isotropic turbulence, show that is has a negligible contribution to 
the evolution of turbulent kinetic energy. Then, an analysis for the case of homogeneous 
turbulence is performed to obtain a model for the non-oscillatory pressure-dilatation. This 
model algebraically relates the pressure-dilatation to quantities traditionally obtained in 
incompressible turbulence closures. The model is validated by direct comparison with the 
pressure-dilatation data obtained from the simulations. 
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Smith, Ralph C., and Kenneth L. Bowers: Sinc-Galerkin estimation of diffusitivity in parabolic 
problems. ICASE Report No. 91-43, May 7, 1991, 36 pages. Submitted to Inverse Problems. 


A fully Sinc-Galerkin method for the numerical recovery of spatially varying diffusion 
coefficients in linear parabolic partial differential equations is presented. Because the pa- 
rameter recovery problems are inherently ill-posed, an output error criterion in conjunction 
with Tikhonov regularization is used to formulate them as infinite-dimensional minimization 
problems. The forward problems are discretized with a sine basis in both the spatial and 
temporal domains thus yielding an approximate solution which displays an exponential con- 
vergence rate and is valid on the infinite time interval. The minimization problems are then 
solved via a quasi-Newton/trust region algorithm. The Z-curve technique for determining 
an appropriate value of the regularization parameter is briefly discussed, and numerical ex- 
amples are given which demonstrate the applicability of the method both for problems with 
noise-free data as well as for those whose data contains white noise. 


Hall, Philip, and Helen Morris: On the instability of boundary layers on heated flat plates. 
ICASE Report No. 91-44 , May 29, 1991, 71 pages. To be submitted to Journal of Fluid 
Mechanics. 


The stability of a boundary layer on a heated flat plate is investigated in the linear 
regime. The flow is shown to be unstable to longitudinal vortex structures which in general 
develop in a nonparallel manner in the streamwise direction. Solutions of the nonparallel 
equations are obtained numerically at 0(1) values of the appropriate stability parameter, 
ie the Grcishof number. The particular cases investigated relate to the situations when the 
instability is induced by localized or distributed wall roughness or nonuniform wall heating. 
The case when the vortices are induced by freestream disturbances is also considered. The 
fastest growing mode is found to be governed by a quasi-parallel theory at high wavenumbers. 

... aven..inber and growth rate of the fastest growing mode are found in closed form. At 
low wavenumbers the vortex instability is shown to be closely related to Tollmein-SchUchting 
wave. , the effect of wall heating or cooling on the latter type of instability is discussed. 


Protzel, Peter W., Daniel L. Palumbo, and Michael K. Arras: Performance and fault- 
tolerance of neural networks for optimization. ICASE Report No. 91-45 , June 6, 1991, 33 
pages. Submitted to IEEE Transactions on Neural Networks. 


One of the key benefits of future hardware implementations of certain Artificial Neural 
Networks (ANNs) is their apparently ’’built-in” fault-tolerance, which makes them potential 
candidates for critical tasks with high reliability requirements. This paper investigates the 
fault-tolerance characteristics of time-continuous, recurrent ANNs that can be used to solve 
optimization problems. The performance of these networks is first illustrated by using well- 
known model problems like the Traveling Salesman Problem and the Assignment Problem. 
The ANNs are then subjected to up to 13 simultaneous ”stuck-at-l” or ”stuck-at-0” faults 
for network sizes of up to 900 ’’neurons.” The effect of these faults on the performance is 
demonstrated and the cause for the observed fault-tolerance is discussed. An application is 
presented in which a network performs a critical ta.sk for a real-time distributed processing 
system by generating new task allocations during the reconfiguration of the system. The 
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performance degradation of the ANN under the presence of fa,ults is investigated by large- 
scale simulations and the potential benefits of delegating a critical task to a fault-tolerant 
network are discussed. 


Bryan, Kurt, and Michael Vogelius: A uniqueness result concerning the identification of a col- 
lection of cracks from finitely many electrostatic boundary measurements. ICASE Report No. 
91-46 . June 11, 1991, 18 pages. Submitted to SIAM Journal of Mathematical Analysis. 

We consider the problem of locating and identifying a collection of finitely many cracks 
inside a planar domain from measurements of the electrostatic boundary potentials induced 
by specified current fluxes. It is shown that a collection of n or fewer cracks can be uniquely 
identified by measuring the boundary potentials induced by n -f 1 specified current fluxes, 
consisting entirely of electrode pairs. 


Banks, H.T., K. Ito, and C. Wang: Exponentially stable approximations of weakly damped 
wave equations. ICASE Report No. 91-47 , June 11, 1991, 40 pages. To appear in Proceed- 
ings of the International Conference on Control of DPS, Voran (July 1990), Birkhauser. 

We consider wave equations with damping in the boundary conditions. Techniques to as- 
certain the uniform preservation under approximation of exponential stability are presented. 
Several schemes for which preservation can be guaranteed are analyzed. Numerical results 
that demonstrate the lack of stability under approximation for several popular schemes (in- 
cluding standard finite difference and finite element schemes) are given. 


Banks, H.T., D. Cioranescu, A. Das, R. Miller and D.A. Rebnord: Homogenization techniques 
and estimation of material parameters in distributed structures. ICASE Report No. 91-48 , 
June 12, 1991, 24 pages. To appear in Proc. 2nd Conf. on Computational Control, 
Birkhauser, (August 1991). 

We discuss the use of homogenization techniques to derive approximate models with sim- 
ple geometry for physical models of grids and trusses which have a complex geometry that 
gives rise to computational difficulties. Our presentation is in the context of inverse or pa- 
rameter estimation problems for composite material structures with unknown characteristics 
such as stiffness and internal damping. We present the necessary theoretical foundations for 
this approach and discuss comparison of modal properties of the resulting homogenization 
model for a two-dimensional grid structure with modal properties observed in experiments 
with this grid. 
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Ou, Yuh-Roung, and John A. Burns: Optimal control of lift/drag ratios on a rotating cylin- 
der. ICASE Report No. 91-49 , June 14, 1991, 13 pages. Submitted to Applied Mathematics 
Letters. 


We present the numerical solution to a problem of maximizing the lift to drag ratio by 
rotating a circular cylinder in a two-dimensional viscous incompressible flow. This problem 
is viewed as a test case for the newly developing theoretical and computational methods 
for control of fluid dynamics systems. We show that the time averaged lift to drag ratio 
for a fixed finite-time interval achieves its maximum value at an optimal rotation rate that 
depends on the time interval. 


Duck, Peter W.: The unsteady laminar boundary layer on an axisymmetric body subject to 
small amplitude fluctuations in the free-stream velocity. ICASE Report No. 91-50 , June 20, 
1991, 43 pages. To appear in J. Fluid Mechanics. 


The effect of small amplitude, time-periodic, freestream disturbances on an otherwise 
steady axisymmetric boundary layer on a circular cylinder is considered. Numerical solutions 
of the problem are presented, and an asymptotic solution, valid far downstream along the 
axis of the cylinder is detailed. Particular emphasis is placed on the unsteady eigensolutions 
that occur far downstream, which turn out to be very different from the analogous planar 
eigensolutions. These axisymmetric eigensolutions are computed numerically and also are 
described by asymptotic analyses valid for low and high frequencies of oscillation. 


Lasseigne, D. G., and T. L. Jackson: Stability of a non- orthogonal stagnation flow to three 
dimensional disturbances. ICASE Report No. 91-51 , June 20, 1991, 28 pages. To appear in 
Theoretical and Computational Fluid Dynamics. 


A similarity solution for a low Mach number non-orthogonal flow impinging on a hot or 
cold plate is presented. For the constant density ca.se, it is known that the stagnation point 
shifts in the direction of the incoming flow and that this shift increases as the angle of attack 
decreases. When the effects of density variations are included, a critical plate temperature 
exists; above this temperature the stagnation point shifts away from the incoming stream as 
the angle is decreased. This flow field is believed to have applications to the reattachment 
zone of certain separated flows or to a lifting body at a high angle of attack. Finally, we 
examine the stability of this non-orthogonal flow to self-similar, three dimensional distur- 
bances. Stability characteristics of the flow are given as a function of the parameters of this 
study: ratio of the plate temperature to that of the outer potential flow and angle of attack. 
In particular, it is shown that the angle of attack can be scaled out by a suitable definition 
of aji equivalent wavenumber and temporal growth rate, and the stability problem for the 
non-orthogonal case is identical to the stability problem for the orthogonal case. By use of 
this scaling, it can be shown that decreasing the angle of attack decreases the wavenumber 
and the magnitude of the temporal decay rate, thus making nonlinear effects important. For 
small wavenumbers it is shown that cooling the plate decreases the temporal decay of the 
least stable mode, while heating the plate has the opposite effect. For moderate to large 
wavenumbers, density variations have little effect except that there exists a range of cool 
plate temperatures for which these disturbances are extremely stable. 
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0v6rina.n, Andrea L., and John Van Rosendale; Mapping implicit spectral methods to dis- 
tributed memory architectures. ICASE Report No. 91-52 , June 28, 1991, 9 pages. To appear 
in Proc. 5th SIAM Conf. on Parallel Computing, Houston, TX, March 1991. 


Spectral methods have proven invaluable in numerical simulation of PDEs, but the fre- 
quent global communication required raises a fundamental barrier to their use on highly 
parallel archtitectures. To explore this issue, we implemented a t^ee dimensional implicR 
spectral method on an Intel hypercube. Utilization of about 50% was achieved on a 32 
node iPSC/860 hypercube, for a 64 x 64 x 64 Fourier-spectral grid; finer grids yield higher 

utilizations. , . • i j u- -j • 

Chebyshev-spectral grids are more problematic, since plane-relaxation based multigrid is 
required. However, by using a semicoarsening multigrid algorithin, and by relaxing all multi- 
grid levels concurrently, relatively high utilizations were also achieved in this harder case. In 
fact, since the amount of work per processor was higher in this case, we achieved somewhat 
higher utilization, typically 60% on moderate sized problems. Thus spectral methods remain 
attractive on the current generation of distributed memory architectures. 


Tanveer, S.: Evolution of a Hele Shaw interface for small surface tension. ICASE Report No. 
91-53 . June 28, 1991, 63 pages. Submitted to Phil. Trans. Roy. Soc. London. 


We present analytical evidence suggesting that for zero surface tension, the initial value 
problem for a Hele-Shaw flow described by a conformal map z{^,t) at time t from inside 
the unit f semi-circle (or unit circle for a radial geometry) to the physical flow domain, 
is well posed when the spatial domain is extended to the unphysical domain [(1 >1 m 
contrast to calculation in the physical domain |C| ^ 1 • particular, we present evidence 
to suggest that there exists solution z{C,t) with branch point singularities behavior of a 
power type, i.e. ~ Ho (C ~ (/^ integral) as > Cs{i) The form of the 

singularity remains invariant in time though its strength Bq changes. A singularity does 
not impinge the physical domain in finite time when ^ > ^ • 

If B denote a nondimensional surface tension, analytical evidence indicates that for 
Q B « 1 for t « surface tension modifies the singularity structure only 

locally for certain kinds of singularities; yet at each initial zero of Zf for initial conditions 
independent of B , new singularities in the higher order terms of the outer asymptotic 
expansion in powers of B are born in 0{B'^'^) time that move towards the physical 
domain. The singularities are smoothed out over an inner scale depending on B ; yet as 
they approach the physical domain, they can correspond to irregular features in a length 
scale larger than a surface tension dependent scale but smaller than a global scale like the 
channel width. 


Dando, Andrew, and Sharon 0. Seddougui: The inviscid compressible Gortler problem. 

ICASE Report No. 91-54 , July 18, 1991, 40 pages. To be submitted to IMA Journal of 

Applied Mathematics. 


In this paper we investigate the growth rates of Gortler vortices in a compressible flow 
in the inviscid limit of large Gortler number. Numerical solutions are obtained for 0(1) 
wavenumbers. The further limits of (i) large Mach number and (ii) large wavenumber with 
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0(1) Mach number are considered. We show that two different types of disturbance modes can 
appear in this problem. The first is a wall layer mode, so named as it has its eigenfunctions 
trapped in a thin layer near the wall. The other mode we investigate is confined to a thin 
layer away from the wall and termed a trapped layer mode for large wavenumbers and an 
adjustment layer mode for large Mach numbers, since then this mode has its eigenfunctions 
concentrated in the temperature adjustment layer. We are able to investigate the near 
crossing of the modes which occurs in each of the limits mentioned. 


Nicol, David M.: Rectilinear partitioning of irregular data parallel computations. ICASE Re- 
port No. 91-55, July 12, 1991, 29 pages. Submitted to Journal of Parallel and Distributed 
Computing. 


This paper describes new mapping algorithms for domain-oriented data-parallel compu- 
tations, where the workload is distributed irregularly throughout the domain, but exhibits 
localized communications patterns. We consider the problem of partitioning the domain for 
parallel processing in such a way that the workload on the most heavily loaded processor is 
minimized, subject to the constraint that the partition be perfectly rectilinear. Rectilinear 
partitions are useful on architectures that have a fast local mesh network and a relatively 
slower global network; these partitions heuristically attempt to maximize the fraction of 
communication carried by the local network. This paper provides an improved algorithm 
for finding the optimal partition in one dimension, new algorithms for partitioning in two 
dimensions, and shows that optimal partitioning in three dimensions is NP-complete. We 
discuss our application of these algorithms to real problems. 


Swanson, R. C., E. Turkel, and J. A. White: An effective multigrid method for high-speed 
flows. ICASE Report No. 91-56 , July 15, 1991, 27 pages. Submitted to Journal of Commu- 
nications in Applied Numerical Methods. 


We consider the use of a multigrid method with central differencing to solve the Navier- 
Stokes equations for high-speed flows. The time-dependent form of the equations is inte- 
grated with a Runge-Kutta scheme accelerated by local time stepping and variable coefficient 
implicit residual smoothing. Of particular importance are the details of the numerical dis- 
sipation formulation, especially the switch between the second and fourth difference terms. 
Solutions are given for the two-dimensional laminar flow over a circular cylinder and a 15 
degree compression ramp. 


Turkel, E., R. C. Swanson, V. N. Vatsa, and J. A. White: Multigrid for hypersonic viscous 
two- and three-dimensional flows. ICASE Report No. 91-57 , July 15, 1991, 23 pages. Pro- 
ceedings of the AIAA 10th Computational Fluid Dynamics Conference, Honolulu, Hawaii, 
June 24-27, 1991, pp. 501-517. 


We consider the use of a multigrid method with central differencing to solve the Navier- 
Stokes equations for hypersonic flows. The time-dependent form of the equations is inte- 
grated with an explicit Runge-Kutta scheme accelerated by local time stepping and implicit 
residual smoothing. Variable coefficients are developed for the implicit process that remove 
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the diffusion limit on the time step, producing significant improvement in convergence. A 
numerical dissipation formulation that provides good shock-capturing capability for hyper- 
sonic flows is presented. This formulation is shown to be a crucial aspect of the multigrid 
method. Solutions are given for two-dimensional viscous flow over a NACA 0012 airfoil and 
three-dimensional viscous flow over a blunt biconic. 


Speziale, Charles G., and Peter S. Bernard: The energy decay in self-preserving isotropic 
turbulence revisited. ICASE Report No. 91-58 , July 17, 1991, 43 pages. To be submitted to 
the Journal of Fluid Mechanics. 


The assumption of self-preservation allows for an analytical determination of the energy 
decay in isotropic turbulence. Batchelor (1948), who was the first to carry out a detailed 
study of this problem, based his analysis on the assumption that the Loitsianskii integral 
is a dynamic invariant — a widely accepted hypothesis that was later discovered to be 
invalid. Nonetheless, it appears that the self-preserving isotropic decay problem has never 
been reinvestigated in depth subsequent to this earlier work. In the present paper^ such an 
analysis is carried out, yielding a much more complete picture of self- preserving isotropic 
turbulence. It is proven rigorously that complete self-preserving isotropic turbulence admits 
two general types of asymptotic solutions: one where the turbulent kinetic energy /C ~ f 
and one where K ~ with an exponent a > 1 that is determined explicitly by the 
initial conditions. By a fixed point analysis and numerical integration of the exact one- 
point equations, it is demonstrated that the A ~ t ^ power law decay is the asymptoticahy 
consistent high-Reynolds-number solution; the K ~ t~°‘ decay law is only achieved i^i. the 
limit as t oo and the turbulence Reynolds nurnber vanishes. Arguments are provided 
which indicate that a A ~ power law decay is the asymptotic state toward which a 
complete self-preserving isotropic turbulence is driven at high Reynolds numbers in order to 
resolve the imbalance between vortex stretching and viscous diffusion. Unlike in previous 
studies, the asymptotic approach to a self-preserving state is investigated which uncovers 
some surprising results. 


Amitai, Dganit, Amir Averbuch, Samuel Itzikowitz, and Eli Turkel: Asynchronous and 
corrected- asynchronous numerical solutions of parabolic PDEs on MIMD multiprocessors. 
ICASE Report No. 91-59 , July 17, 1991, 23 pages. Submitted to SIAM Journal on Scientific 

and Statistical Computing. 


A major problem in achieving significant speed-up on parallel machines is the overhead 
involved with synchronizing the concurrent processes. Removing the synchronization con- 
straint has the potential of speeding up the computation. We present asynchronous (AS) 
and corrected-asynchronous (CA) finite difference schemes for the multi-dimensional heat 
equation. Although our discussion concentrates on the Euler scheme for the solution of the 
heat equation, it has the potential of being extended to other schemes and other parabolic 
PDEs. These schemes are analyzed and implemented on the shared-memory multi-user 
Sequent Balance machine. Numerical results for one and two dimensional problems are pre- 
sented. It is shown experimentally that synchronization penalty can be about 50% of run 
time: in most cases, the asynchronous scheme runs twice as fast as the parallel synchronous 
scheme. In general, the efficiency of the parallel schemes increases with processor load, with 
the time-level, and with the problem dimension. The efficiency of the AS may reach 90% 
and over, but it provides accurate results only for steady-state values. The CA, on the other 
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hand, is less efficient but provides more accurate results for intermediate (non steady-state) 
values. 


Jacobs, P.A.: Simulation of transient flow in a shock tunnel and a high Mach number noz- 
zle. ICASE Report No. 91-60 , July 18, 1991, 15 pages. To be presented at the 4th Int. 
Symposium on Computational Fluid Dynamics, Davis, CA, September 1991. 


A finite- volume Navier-Stokes code was used to simulate the shock-reflection and nozzle 
starting processes in an axisymmetric shock tube and a high Mach number nozzle. The 
simulated nozzle starting processes were found to match the classical quasi-one-dimensional 
theory and some features of the experimental measurements. The shock-reflection simulation 
illustrated a new mechanism for the driver-gas contamination of the stagnated test gas. 


Thangam, S.: Analysis of two-equation turbulence models for recirculating flows. ICASE Re- 
port No. 91-61, July 19, 1991, 22 pages. Submitted to Journal of Fluids Engineering. 


The two-equation K — e model is used to analyze turbulent separated flow past a 
backward-facing step. It is shown that if the model constants are modified to be consis- 
tent with the accepted energy decay rate for isotropic turbulence, the dominant features of 
the flow field - namely - the size of the separation bubble and the streamwise component 
of the mean velocity, can be accurately predicted. In addition, except in the vicinity of the 
step, very good predictions for the turbulent shear stress, the wall pressure and the wall 
shear stress are obtained. The model is also shown to provide good predictions for the tur- 
bulence intensity in the region downstream of the reattachment point. Estimated long-time 
growth rates for the turbulent kinetic energy and dissipation rate of homogeneous shear flow 
are utilized to develop an optimal set of constants for the two equation K — € model. The 
physical implications of the model performance are also discussed. 


Nessyahu, Haim, and Eitan Tadmor: The convergence rate of approximate solutions for 
nonlinear scalar conservation laws. ICASE Report No. 91-62 , July 24, 1991, 21 pages. Sub- 
mitted to SIAM Journal on Numerical Analysis. 


Let {u'(x, ^)}£>o be a family of approximate solutions for the nonlinear scalar conservation 
law Ut -\- fx{u) = 0 with C^initial data. Assume that {u'(a;,i)} are Xip+-staMe in the sense 
that they satisfy Oleinik’s E-entropy condition. We prove that if these approximate solutions 
are Ltp'- consistent , i.e., if l|u'(-,0) - u(-, 0)l|i,ip/(i) -|- ||u' + fx{v‘)\\Lip>(x,t) = then they 

converge to the entropy solution and the convergence rate estimate, ||u'(-, 0 = 

0(e), holds. Consequently, sharp L’’ and pointwise error estimates are derived. 

We demonstrate these convergence rate results in the context of entropy satisfying finite- 
difference and Glimm’s schemes. 
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Seddougui, Sharon 0., and Andrew P. Bassom: Nonlinedv instability of hypersonic flow past 
a wedge, ICASE Report No. 91-63 , July 25, 1991, 42 pages. To be submitted to Journal of 
Theoretical and Computational Fluid Dynamics. 


The nonlinear stability of a compressible flow past a wedge is investigated in the hyper- 
sonic limit. The analysis follows the ideas of a weakly nonlinear approach as first detailed 
by Smith (1979). Interest is focussed on Tollmien-Schlichting waves governed by a triple 
deck structure and it is found that the attached shock can profoundly affect the stability 
characteristics of the flow. In particular, it is shown that nonlinearity tends to have a sta- 
bilising influence. The nonlinear evolution of the Tollmien-Schlichting mode is described in 
a number of asymptotic limits which were first identified by Cowley and Hall (1990) in their 
linearised account of the current problem. 


Duck, Peter W., Gordon Erlebacher, and M. Yousuff Hussaini: On the linear stability of com- 
pressible plane Couette flow, ICASE Report No. 9T64 , July 29, 1991, 52 pages. Submitted 

to Journal of Fluid Mechanics. 


The linear stability of compressible plane Couette flow is investigated. The correct and 
proper basic velocity and temperature distributions are perturbed by a small amplitude nor- 
mal mode disturbance. The full small amplitude disturbance equations are solved numeri- 
cally at finite Reynolds numbers, and the inviscid limit of these equations is then investigated 
in some detail. It is found that instability can occur, although the stability characteristics of 
the flow are quite different from unbounded flows. The effects of viscosity are also calculated, 
asymptotically, and shown to have a stabilizing role in all the cases investigated. Exceptional 
regimes to the problem occur when the wavespeed of the disturbances approaches the ve- 
locity of either of the walls, and these regimes are also analyzed in some detail. Finally, the 
effect of imposing radiation-type boundary conditions on the upper (moving) wall (in place 
of impermeability) is investigated, and shown to yield results common to both bounded and 
unbounded flows. 


Yakhot, V., S. Thangam, T.B. Gatski, S.A. Orszag, and C.G. Speziale: Development of 
turbulence models for shear flows by a double expansion technique, ICASE Report No, 91-65 , 
July 29, 1991, 26 pages. Submitted to Physics of Fluids A. 


Turbulence models are developed by supplementing the renormalization group (RNG) 
approach of Yakhot & Orszag with scale expansions for the Reynol^ ^ress and production 
of dissipation terms. The additional expansion parameter {rj = SKfS) is the ratio of the 
turbulent to mean strain time scale. While low-order expansions appear to provide an 
adequate description for the Reynolds stress, no finite truncation of the expansion for the 
production of dissipation term in powers of tj suffices — terms of all orders must be retained. 
Based on these ideas, a new two-equation model and Reynolds stress transport model are 
developed for turbulent shear flows. The models are tested for homogeneous shear flow 
and flow over a backward facing step. Comparisons between the model predictions and 
experimental data are excellent. 
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Vatsa, Veer N., Eli Turkel, and J. S. Abolhassani: Extension of multigrid methodology to su- 
personic/hypersonic 3-D viscous flows, ICASE Report No. 91-66 , August 2, 1991, 24 pages. 
Submitted to Internation Journal of Numerical Methods in Fluid Mechanics. 


A multigrid acceleration technique developed for solving the three-dimensional Navier- 
Stokes equations for subsonic/transonic flows has been extended to supersonic/hypersonic 
flows. An explicit multistage Runge-Kutta type of time-stepping scheme is used as the basic 
algorithm in conjunction with the multigrid scheme. Solutions have been obtained for a 
blunt conical frustum at Mach 6 to demonstrate the applicability of the multigrid scheme 
to high-speed flows. Computations have also been performed for a generic High-Speed Civil 
Transport configuration designed to cruise at Mach 3, These solutions demonstrate both the 
efficiency and accuracy of the present scheme for computing high-speed viscous flows over 
configurations of practical interest. 


Ou, Yuh-Roung: Control of oscillatory forces on a circular cylinder by rotation. ICASE Re- 
port No. 91-67, August 15, 1991, 17 pages. To appear in Proceedings of 4th International 
Symposium on Computational Fluid Dynamics... Davis, CA, September 1991. 


The temporal development of forces acting on a rotating cylinder is investigated numer- 
ically in response to a variety of time-dependent rotation rates. Solutions are presented for 
several types of rotation that illustrate significant effects of the rotation rate on lift, drag and 
lift/drag coefficients. Of special interest is the formulation of an optimal control problem for 
the case of constant speed of rotation. We find an optimal rotation rate that achieves the 
maximum value of time-averaged list /drag ratio. 


Geer, James F., and Carl M. Andersen: Resonant frequency calculations using a hybrid 
perturbation- Galerkin technique. ICASE Report No. 91-68 , September 11, 1991, 36 pages. 
To appear in Journal of Applied Mechanics. 


A two step hybrid perturbation-Galerkin technique is applied to the problem of deter- 
mining the resonant frequencies of one- or several-degree(s)-of-freedom nonlinear systems 
involving a parameter. In step one, the Lindstedt-Poincare method is used to determine 
perturbation solutions which are formally valid about one or more special values of the 
parameter (e.g. for small or large values of the parameter). In step two, a subset of the 
perturbation coordinate functions determined in step one is used in a Galerkin type approx- 
imation. The technique is illustrated for several one-degree-of-freedom systems, including 
the Duffing and van der Pol oscillators, as well as for the compound pendulum. For all of the 
examples considered, it is shown that the frequencies obtained by the hybrid technique using 
only a few terms from the perturbation solutions are significantly more accurate than the 
perturbation results on which they are based, and they compare very well with frequencies 
obtained by purely numerical methods. 
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Macaraeg, Michele G., T. L. Jackson and M. Y. Hussaini: Ignition and structure of a laminar 
diffusion flame in the field of a vortex. ICASE Report No. 91-69 , September 30, 1991, 35 
pages. Submitted to Combustion Science of Technology. 


The distortion of flames in flows with vortical motion is examined via asymptotic analysis 
and numerical simulation. The model consists of a constant-density, one-step, irreversible 
Arrhenius reaction between initially unmixed species occupying adjacent half-planes which 
are then allowed to mix and react in the presence of a vortex. The evolution in time of the 
temperature and mass-fraction fields is followed. Emphasis is placed on the ignition time 
and location as a function of vortex Reynolds number and initial temperature differences 
of the reacting species. The study brings out the influence of the vortex on the chemical 
reaction. In all phases, good agreement is observed between asymptotic analysis and the full 
numerical solution of the model equations. 


Naik, Naomi H., and John Van Rosendale: The improved robustness of multigrid elliptic 
solvers based on multiple semicoarsened grids. ICASE Report No. 91-70 , September 6, 1991, 
18 pages. Submitted to SIAM Journal of Numerical Analysis. 


Multigrid convergence rates degenerate on problems with stretched grids or anisotropic 
operators, unless one uses line or plane relaxation. For three dimensional problems, only 
plane relaxation suffices, in general. While line and plane relaxation algorithms are efficient 
on sequential machines, they are quite awkward and inefficient on parallel machines. This 
paper presents a new multigrid algorithm, based on the use of multiple coarse grids, that 
eliminates the need for line or plane relaxation in anisotropic problems. We develop this 
algorithm, and extend the standard multigrid theory to establish rapid convergence for this 
class of algorithms. The new algorithm uses only point relaxation, allowing easy and efficient 
parallel implementation, yet achieves robustness and convergence rates comparable to line 
and plane relaxation multigrid algorithms. 

The algorithm described here is a variant of Mulder’s multigrid algorithm [5] for hyper- 
bolic problems. The latter uses multiple coarse grids to achieve robustness, but is unsuitable 
for elliptic problems, since its V-cycle convergence rate goes to one as the number of levels 
increases. The new algorithm combines the contributions from the multiple coarse grids via 
a local “switch,” based on the strength of the discrete operator in each coordinate direction. 
This improvement allows^us to show that the V-cycle convergence rate is uniformly bounded 
away from one, on model anisotropic problems. Moreover, the new algorithm can be com- 
bined with the idea of concurrent iteration on all multigrid levels to yield a highly parallel 
algorithm for strongly anisotropic problems. 


Carpenter, Mark H., David Gottlieb, and Saul Abarbanel: The stability of numerical bound- 
ary treatments for compact high-order finite-difference schemes. ICASE Report No. 91-71 , 
September 10, 1991, 57 pages. To be submitted to Journal of Computational Physics. 


The stability characteristics of various compact fourth- and sixth-order spatial opera- 
tors are assessed using the theory of Gustafsson, Kreiss and Sundstrom (G-K-S) for the 
semi-discrete Initial Boundary Value Problem (IBVP). These results are then generalized to 
the fully discrete case using a recently developed theory of Kreiss. In all cases, favorable 
comparisons are obtained between G-K-S theory, eigenvalue determination, and numerical 
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simulation. The conventional definition of stability is then sharpened to include only those 
spatial discretizations that are asymptotically stable (bounded, Left Half-Plane eigenval- 
ues). It is shown that many of the higher-order schemes which are G-K-S stable are not 
asymptotically stable. A series of compact fourth- and sixth-order schemes, which are both 
asymptotically and G-K-S stable for the scalar case, are then developed. 


Chapman, Barbara, Piyush Mehrotra, and Hans Zima; Vienna fortran - A fortran language 
extension for distributed memory multiprocessors. ICASE Report No. 91-72 , September 10, 
1991, 29 pages. To appear in book “Languages, Compilers and Runtime Environments for 
Distributed Memory Machines”, Edited by J. Saltz and P. Mehrotra, Elsevier Press. 


Exploiting the performance potential of distributed memory machines requires a careful 
distribution of data across the processors. Vienna Fortran is a language extension of Fortran 
which provides the user with a wide range of facilities for such mapping of data structures. 
However, programs in Vienna Fortran are written using global data references. Thus, the user 
heis the advantages of a shared memory programming paradigm while explicitly controlling 
the placement of data. In this paper, we present the basic features of Vienna Fortran along 
with a set of examples illustrating the use of these features. 


Das, Raja, Ravi Ponnusami, Joel Saltz and Dimitri Mavriplis: Distributed memory compiler 
methods for irregular problems - Data copying reuse and runtime partitioning. ICASE Report 
No. 91-73 . September 17, 1991, 35 pages. To appear in book “Languages, Compilers and 
Runtime Environments for Distributed Memory Machines”, Editors: J. Saltz and P. Mehro- 
tra, Elsevier Press. 


This paper outlines two methods which we believe will play an important role in any 
distributed memory compiler able to handle sparse and unstructured problems. We describe 
how to link runtime partitioners to distributed memory compilers. In our scheme, pro- 
grammers can implicitly specify how data and loop iterations are to be distributed between 
processors. This insulates users from having to deal explicitly with potentially complex 
algorithms that carry out work and data partitioning. 

We also describe a viable mechanism for tracking and reusing copies of off-processor 
data. In many programs, several loops access the same off-processor memory locations. As 
long as it can be verified that the values cissigned to off-processor memory locations remain 
unmodified, we show that we can effectively reuse stored oflF-processor data. We present 
experimental data from a 3-D unstructured Euler solver run on an iPSC/860 to demonstrate 
the usefulness of our methods. 


Otto, S.R.: Stability of the flow around a cylinder: The spin-up problem. ICASE Report No. 

91-74 . September 17, 1991, 25 pages. Submitted to Theoretical and Computational Fluid 
Dynamics. 


Our concern is with the flow around an infinite cylinder, which at a certain instant is 
impulsively staxted to spin. The growth of vortices in the resulting boundary layer occurring 
outside the cylinder is investigated. This layer is essentially a Rayleigh layer which grows 
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with time, so the mechanism involved is similar to that studied in Hall (198^. Vortices 
with wavenumber comparable to the layer thickness are shown to be described by partial 
differential equations. It is found necessary to solve the unsteady partial differential equaticms 
that govern the system numerically. We assume that the Rayleigh layer is thin so particles 
are confined to move in a path with radius of curvature the same as the cylinder. The Gortler 
number is a function of time, so we consider the time scale which produces an order one 
Gortler number. We consider the right hand branch calculation by letting the time tend to 
infinity, also inviscid Gortler modes are considered. 


Jacobs, P.A.: An approximate Riemann solver for hypervelocity flows. ICASE Report No. 
91-75 . September 23, 1991, 15 pages. Submitted as “Technical Notes” to the AIAA Journal. 


We describe an approximate Riemann solver for the computation of hypervelocity flows 
in which there are strong shocks and viscous interactions. The scherne has three stages, 
the first of which computes the intermediate states assuming isentropic waves. A second 
stage, based on the strong shock relations, may then be invoked if the pressure jump across 
either wave is large. The third stage interpolates the interface state from the two initial 
states and the intermediate states. The solver is used as part of a finite- volunie code and is 
demonstrated on two test cases. The first is a high Mach number flow over a sphere while 
the second is a flow over a slender cone with an adiabatic boundary layer. In both cases the 
solver performs well. 


Harten Ami, and Sukumar R. Chakravarthy: Multi-dimensional ENO schemes for general 
geometries. ICASE Report No. 91-76 , September 20, 1991, 70 pages. To be submitted to 

JCP. 


In this paper we present a class of EN 0 schemes for the numerical solution of multidimen- 
sional hyperbolic systems of conservation laws in structured and unstructured grids. This 
is a clciss of shock-capturing schemes which are designed to compute cell-averages to high- 
order of accuracy. The ENO scheme is composed of a piecewise-polynomial reconstruction 
of the solution from its given cell- averages, approximate evolution of the resulting initial- 
value problem, and averaging of this approximate solution over each cell. The reconstruction 
algorithm is based on an adaptive selection of stencil for each cell so as to avoid spurious 
oscillations near discontinuities while achieving high order of accuracy away from them. 


Harten, Ami: Multi-resolution analysis for ENO schemes. ICASE Report No. 91-77 , Septem- 
ber 26, 1991, 16 pages. To appear in Proceedings of the ICASE/LaRC workshop on Algo- 
rithmic Trends in CFD for the 90’s. 


Given a function u{x) which is represented by its cell-averages in cells which are formed 
by some unstructured grid, we show how to decompose the function into various scales of 
variation. This is done by considering a set of nested grids in which the given grid is the 
finest, and identifying in each locality the coarsest grid in the set from which u{x) can be 
recovered to a prescribed accuracy. 
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We apply this multi-resolution analysis to ENO schemes in order to reduce the number of 
numerical flux computations which is needed in order to advance the solution by one time- 
step. This is accomplished by decomposing the numerical solution at the beginning of each 
time-step into levels of resolution, and performing the computation in each locality at the 
appropriate coarser grid. We present an efficient algorithm for implementing this program 
in the one-dimensional case; this algorithm can be extended to the multi-dimensional case 
with cartesian grids. 
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ICASE INTERIM REPORTS 


Das, Subhendu, Joel Saltz, and Harry Berryman: A manual for PARTI runtime primitives 
- Revision 1., Interim Report No. 17 , April 26, 1991, 53 pages. 


Primitives are presented that are designed to help users efficiently prc^ram irregular 
problems (e.g. unstructured mesh sweeps, sparse matrix codes, adaptive mesh partial dmer- 
ential equations solvers) on distributed memory machines. These primitives are also designed 
for use in compilers for distributed memory multiprocessors. Communications patterns are 
captured at runtime, and the appropriate send and receive messages are automatically gen- 

erated. 


Jacobs, P. A.; Single-block Navier-Stokes integrator.. Interim Report No. 18 , July 12, 1991, 
71 pages. 


This report describes a program for the time-integration of the Navier Stokes equations 
on a two-dimensional structured mesh. The flow geometry may be eijer planar or axisym- 
metric. The unusual features of this program are that it is writ^n in C and makes extensive 
use of sophisticated data structures to encapsulate the data. The idea of writing the code 
this way is to make it easier Hhan traditional FORTRAN codes) to parallelize for the 

Multiple-Instruction-Multiple-Data style of parallel computer. 

The integral form of the governing equations are given for cartesian coordinates and then 
the particular discretization used in the code is described. A derivation of the axisymmetric 
equations is given in an appendix. The full version of the code describes a flow domain as a set 
of abutting blocks, each consisting of a tensor-product mesh of quadrilateral cells. However, 
this report considers only the single-block version of the code. The flow field is recorded 
as cell-average values at cell centres and explicit time stepping is used to update conserved 
quantities. MUSCL-type interpolation and a three-stage Riemann solver are used to calculate 
mviscid fluxes across cell faces while central differences (via the divergence theorem) are used 
to calculate the viscous fluxes. The Riemann solver is suitable for flows with very strong 
shocks and does not require the entropy fix as applied to the Roe-type solvers Because 
the code is intended to be a test-bed for implementation on parallel computers, the coding 

details are described in some detail. _ ^ . i j j 

A set of test problems is also included. These exercise various parts of the code and 
should be useful for both validation and performance measurements of the (future) parallel 
implementations. 


Perkins, A. Louise, and Jeffrey S. Scroggs: Proceedings for the ICASE workshop on hetero- 
geneous boundary conditions., Interim Report No. 19 , August 15, 1991, 78 pages. 


Domain Decomposition is a complex problem with many interesting aspects. The choice 
of decomposition can be made based on many different criteria, and the choice of interface 
of internal boundary conditions are numerous. Even more interesting froin a modeling 
perspective is that the various regions under study may have different dynamical balances, 
indicating that different physical processes are dominating the flow in those regions. It may 
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be desirable to use different numerical approximations in the regions where the physical 
processes are dominated by different balances. 

The Institute for Computer Applications in Science and Engineering (ICASE), recogniz- 
ing the need to more clearly define the nature of these complex problems, sponsored this 
workshop on Heterogeneous Boundary Conditions at the NASA Langley Research Center 
in Hampton, Virginia. This proceedings is an informal collection of the presentations and 
discussion groups. It also includes a bibliography that contains many of the references that 
discuss related topics. 

The proceedings begins with summaries of the discussion groups. Then papers describing 
the talks are presented. Lastly, the bibliography is included, and an index by subject is 
provided. 
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ICASE COLLOQUIA 

April 1, 1991 - September 30, 1991 

Name/Affiliation/Title Date 

Padma Raghavan, The Pennsylvania State University April 8 

“Distributed Sparse Matrix Factorization: 

QR and Cholesky Decompositions” 

Marten T. Landahl, Massachusetts Institute of Technology April 9 

“Algebraic Instability and the Near- Wall Reynolds Stresses” 

Hans P. Zima, University of Vienna, Austria April 15 

“Compiling for Distributed-Memory Machines: 

Experiences and Recent Research” 

Ricardo D. Pantazis, Duke University April 16 

“Parallel Solution of the Symmetric Generalized 
Eigenvalue Problem” 

Francois Irigoin, Centre de Recherche en Informatique, France April 19 
“Automatic Interprocedural Parallelization” 

Suresh Chittor, Michigan State University April 19 

“Communication Performance of Multicomputers” 

Alan Sussman, Carnegie-Mellon University May 2 

“Model-Driven Mapping onto a Distributed Memory 
Parallel Computer” 

Michael Gerndt, University of Vienna May 3 

“Compiling Sequential Programs for Shared-Memory and 
Distributed-Memory Multiprocessors” 

Jacques Periaux, Dassault Aviation, France May 9 

“Finite Element Simulations of Hypersonic Reacting Flows: 

Numerical Tools for the Design of the European Hermes 
Space Vehicle” 

Hermann F. Fasel, The University of Arizona, Tucson May 15 

“Numerical Simulation of Laminar- Turbulent Transition 
in Shear Flows” 
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Name/AfRliation/Title 


Date 


Steven Zenith, Ecole Nationale Superieure des Mines de Paris, France 
“Process Interaction Models: A Brief Description of 
and Rationale for Ease” 

Matthew Rosing, University of Colorado, Boulder 
“Language Support for Numerical Programs 
on Distributed Memory Multiprocessors” 

William T. Ashurst, Sandia National Laboratories 
“Geometry of Premixed Flames in Three-Dimensional Turbulence” 

Fabio Bertolotti, Princeton University 

“Investigations of Laminar-Turbulent Transition in Subsonic 
and Supersonic Boundary Layers Using the Parabolic 
Stability Equations” 

Asok Ray, The Pennsylvania State University 

“Damage-Mitigating Control of Aerospace Structures for High 
Performance and Extended Life” 

Rajesh Aggarwal, Columbia University 

“Parallel Algorithms in Computational Mechanics” 

J. Peraire, Imperial College of Science and Medicine, London 
“The Computation of Three Dimensional Flows Using 
Unstructured Meshes” 

Carl Scarbnick, Scientific Computing Associates 
“Algorithms for Solving Problems in Scientific Computation” 

Lawrence Crowl, University of Rochester 
“Architectural Adaptability in Parallel Programming” 

Jean Paul Boujot, CISI Ingenierie, France 
“Will Computer Architectures be Specified by Numerical Methods?” 

I. G. Angus, Northrop Research and Technology Center 
“Parallelism, Object Oriented Programming Methods, 

Portable Software and C-1-+ (A CFD Example)” 


May 16 

May 17 

May 17 
May 22 

June 4 

June 14 
June 19 

June 24 
June 28 
July 15 
July 16 
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Name/AfRliation/Title 


Date 


Mark Jones, Argonne National Laboratory 

“An Efficient Parallel Iterative Solver for Large Sparse Linear Systems” 

Ishfaq Ahmad, Syracuse University 
“Dynamic Task Scheduling and Load Balancing 
for Large Multicomputer Systems” 

Isaac Harari, Stanford University 
“Finite Element Methods for the Reduced Wave Equation: 

Design and Analysis” 

Robert Krasny, University of Michigan 

“Vortex Sheet Computations: Roll-Up, Wakes, Separation” 

John Crepeau, University of Utah 

“Dynamical Systems Study of Mixing and Boundary Layers” 

Daniel Joseph, University of Minnesota 

“Finite Size and Nonlinear Effects in Fluidized Suspensions” 

M. Karpel, Technion - Israel Institute of Technology, Haifa, Israel 
“Multidisciplinary Optimization of Aeroservoelastic Systems 
Using Reduced-Size Models” 

James Antaki, University of Pittsburgh 

“The Artificial Heart: Prospective Need for CFD” 

Charles Fineman, NASA Ames Research Center 
“AIMS: A Software Visualization Tool” 

Philip Rosenau, Technion, Haifa, Israel 
“On Mathematical Modeling of High Gradient Phenomena” 

Rolf Radespiel, DLR, Braunschweig, Germany 
“Progress with Multigrid Schemes for Hypersonic Flow Problems” 

Robert Weaver, University of Colorado 
“Supporting Dynamic Data Structures at the Language Level 
on Distributed Memory Machines” 


July 18 
July 19 

July 22 

August 1 
August 9 
August 13 
August 15 

August 26 
August 28 
August 29 
September 5 
September 10 
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Name/ Affiliation/Title 


Date 


V. Venkatakrishnan, Computer Sciences Corporation 

and ^JASA Ames Research Center September 19 

“A MIMD Implementation of a Parallel Euler Solver 
for Unstructured Grids” 

Ami Harten, Tel Aviv University and ICASE September 20 

“Multi-Dimensional Essentially Non-Oscillatory (ENO) Schemes 
for General Geometries” 

David Mount, University of Maryland September 20 

“Computational Geometry and Mesh Generation” 


51 


ICASE SUMMER ACTIVITIES 


The summer program for 1991 included the following visitors: 


NAME/AFFILIATION 

DATE OF VISIT 

AREA OF INTEREST 

Saul S. Abarbanel 

Tel- Aviv University, Israel 

6/24 - 8/30 

Computational Fluid Dynamics 

Remi Abgrall 
INRIA 

7/1 - 8/30 

Computational Fluid Dynamics on 
Unstructured Meshes 

H. Thomas Banks 

University of Southern California 

6/3 - 6/14 

Control Theory 

Alvin Bayliss 
Northwestern University 

6/10 - 6/13 

Numerical Solution of PDF’s 

Scott Beryman 
Yale University 

5/20 - 8/30 

Software for Parallel Systems 

Shahid Bokhari 

Pakistan University of Engineering 
and Technology 

7/8 - 9/6 

Parallel Computing Systems 

John A. Burns 

VPI and State University 

6/10-14 

Control Theory 

Claudio Canuto 
Politecnico di Torino, Italy 

9/10-20 

Spectral Methods for PDE’s 

Barbara Chapman 
University of Vienna 

8/14-28 

Compiler Development 
for Multiprocessors 

Craig Chase 
Cornell University 

5/20 - 8/30 

Distributed Memory Architectures 

William Criminale 
University of Washington 

8/26 - 9/6 

Transition and Turbulence 

Wai-Sun Don 
Brown University 

6/10 - 7/26 

Computational Fluid Dynamics 
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NAME/AFFILIATION 


DATE OF VISIT AREA OF INTEREST 


Peter L. Duck 

University of Manchester, England 


7/8 - 9/20 Computational Fluid Dynamics 


Dennis Gannon 8/5-9 Parallel Programming Environments 

Indiana University, Bloomington 


James Geer 

State University of New York 


9/19-23 Asymptotic Methods 

for Differential Equations 


Ann Gelb 
Brown University 


6/17-28 Numerical Analysis 


David Gottlieb 6/17 - 8/16 Numerical Methods for PDE’s 

Brown University 


Satyanarayan Gupta 
Old Dominion University 


6/3 - 8/12 Unstructured Mesh Calculations 
for the Connection Machine 


Philip Hall 6/24 - 8/2 Computational Fluid Dynamics 

University of Manchester, England 9/25-10/18 


Amiram Harten 9/1-27 Numerical Methods for PDE’s 

Tel-Aviv University 


Daniel J. Inman 6/3-14 

State University of New York, Buffalo 


Numerical Methods for 
Inverse Problems in Distributed 
Systems 


Kazufumi Ito 6/3-7 Control Theory 

University of California, Los Angeles 


Ashwajii Kapila 

Rensselaer Polytechnic Institute 


8/5-16 Mathematical Aspects 

of Combustion Processes 


Mordechay Karpel 
Technion-IIT, Haifa, Israel 


7/8 - 9/6 


Modeling of Aeroservoelastic 
Systems 


David R. Keyes 
Yale University 


8/12-30 Parallel Numerical Procedures 
Combustion 


Fumio Kojima 

Osaka Institute of Technology, Japan 


8/5-30 Probabilistic and Stochastic Methods 

for Optimal Control Problems 
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NAME/AFFILIATION 

DATE OF VISIT 

AREA OF INTEREST 

Robert Krasny 
University of Michigan 

7/29 - 8/2 

Vortex Methods and Vortex 
Dynamics 

Marten Landahl 

Massachusetts Institute of Technology 

8/26-30 

Transition and Turbulence 

Jacques Liandrat 
I.M.S.T., France 

6/10-14 

Numerical Algorithms for Nonlinear 
PDE 

Robert MacCormack 
Stanford University 

7/29 - 8/2 

Numerical Methods for the Solution 
of the Equations of Fluid Mechanics 

Moshe Matalon 
Northwestern University 

6/3-14 

Hydrodynamic Instability 
in Chemically Reacting 
Flows 


i 


I 


K. W. Morton 
Oxford University 

9/13-21 

H. S. Mukunda 

Indian Institute of Science, Bangalore 

5/29 - 6/21 

Naomi H. Naik 
Vassar College 

8/5-30 

David Nicol 

College of William and Mary 

5/13-17 

8/16-22 

D. Papageorgiou 

New Jersey Institute of Technology 

5/20-31 

Terry Pratt 
University of Virginia 

5/21-24 

6/11-14 

6/18-21 

8/5-9 

8/26-30 

Rolf Radespiel 

DLR, Institute for Design Aerodynamics 

6/17 - 9/14 

Asok Ray 

5/20 - 6/7 


The Pennsylvania State University 


Computational Fluid Dynamics 
High Speed Reacting Flows 
Multigrid Methods 


Techniques for Mapping Algorithms 
onto Parallel Systems 

Computational Fluid Dynamics 


Characteristics of Languages 
for Parallel Computers 


Numerical Methods for 
High Speed Flows 

Control Systems for Flexible 
Space Structures 
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NAME/AFFILIATIQN 


DATE OF VISIT AREA OF INTEREST 


Philip Roe 

The University of Michigan 

8/5-23 

Computational Fluid Dynamics 

Jeffrey Scroggs 

North Carolina State University 

5/10 & 17 
8/5-16 

Domain Decomposition Techniques 
for Partial Differential 
Equations 

Chi- Wang Shu 
Brown University 

6/3 - 7/5 

Computational Fluid Dynamics 

Alex Solomonoff 
Brown University 

6/13-28 

Numerical Methods for PDF’s 

Eitan Tadmor 
Tel-Aviv University 

7/8 - 8/2 

Numerical Solution of 
Partial Differential Equations 

Saleh Tanveer 
Ohio State University 

6/17-28 

9/9-13 

Modelling of Crystal Growth 

Hein Tran 

North Carolina State University 

6/3-7 

Control Theory 

Eli Turkel 
Tel-Aviv University 

6/24 - 8/23 

Computational Fluid Dynamics 

Bram van Leer 

University of Michigan, Ann Arbor 

8/12-13 

9/16-20 

Computational Fluid Dynamics 

Sesh Venugopal 
Rutgers University 

5/20 - 8/16 

Partitioning of Scientific 
Problems on Parallel Computers 

Eric Voth 
Brown University 

6/24 - 7/8 

Computational Fluid Dynamics 

Claus Wagner 
University of Florida 

6/10 - 7/6 

Problems of Instabilities 
in Incompressible Flows 

Yun Wang 

University of Southern California 

6/10-14 

Mini Ma^t Models for Identification 
and Control 

Hans Zima 
University of Vienna 

8/14-28 

Compiler Development for 
Multiprocessors 
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OTHER ACTIVITIES 


The Transition and Turbulence Workshop co-sponsored by ICASE and NASA Langley 
Research Center was held July 8 - August 2,1991, at NASA Langley Research Center. Over 
one hundred people attended this workshop. The objective was to encourage development 
and validation of Reynolds-averaged two-equation and second-order closure models as well 
as subgrid scale models for the transitional and fully turbulent zone. A volume of the 
proceedings from this conference will be published by Springer- Verlag in the near future. 


On September 15-17, 1991 a Workshop on Algorithmic Trends in CFD for the 90’s co- 
sponsored by ICASE and NASA Langley Research Center was held at the Holiday Inn- 
Hampton- Coliseum Hotel and Conference Center. One hundred and ten people attended 
this workshop. The purpose of the workshop was to bring together numerical analysts and 
computational fluid dynamicists i) to assess the state of the art in the areas of numerical 
analysis particularly relevant to CFD, ii) to identify promising new developments in various 
areas of numerical analysis that will have impact on CFD, and iii) to establish a long-term 
prospective focussing on opportunities and needs. 

A volume of the proceedings from this workshop will be published by Springer- Verlag in 

the near future. 
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ICASE STAFF 


I. ADMINISTRATIVE 

Robert G. Voigt, Director Ph.D., Mathematics, University of Maryland, 1969. Numerical 
Algorithms for Parallel Computers 

Linda T. Johnson, Office and Financial Administrator 

Judy L. Batten, Office Assistant (Part-time) 

Etta M. Blair, Personnel/Bookkeeping Secretary 

Barbara A. Cardasis, Administrative Secretary 

Holly D. Joplin, Office Assistant (Part-time) 

Rosa H. Milby, Short-term Housing/Office Secretary 

Shelly D. Millen, Technical Publications Secretary 

Barbara R. Stewart, Office Assistant (Part-time) 

Emily N. Todd, Executive Secret ary /Visitor Coordinator 

IL SCIENCE COUNCIL for APPLIED MATHEMATICS and COMPUTER 
SCIENCE 

Tony Chan, Professor, Department of Mathematics, University of California at Los Angeles, 
S. Lennart Johnsson, Thinking Machines Corporation. 

Michael O’Donnell, Professor, Department of Computer Science, University of Chicago. 

Joseph Oliger, Professor, Computer Science Department, Stanford University. 

Robert O’Malley, Jr., Chairman, Department of Mathematical Sciences, Rensselaer Poly- 
technic Institute. 
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Stanley J. Osher, Professor, Mathematics Department, University of California. 

John Rice, Chairman, Department of Computer Science, Purdue University. 

Burton Smith, Tera Computer Company, Seattle, WA. 

Robert G. Voigt, Director, Institute for Computer Applications in Science and Engineering, 
NASA Langley Research Center. 

III. ASSOCIATE MEMBERS 

Saul S. Abarbanel, Professor, Department of Applied Mathematics, Tel-Aviv University, 
Israel. 

H. Thomas Banks, Professor, Center for Applied Mathematical Sciences, University of South- 
ern California. 

David Gottlieb, Professor, Division of Applied Mathematics, Brown University. 

Peter D. Lax, Professor, Courant Institute of Mathematical Sciences, New York University. 
Merrell L. Patrick, New Technologies Program Director, National Science Foundation. 

IV. CHIEF SCIENTIST 

M. Yousuff Hussaini - Ph.D., Mechanical Engineering, University of California, 1970. Com- 
putational Fluid Dynamics. (Beginning April 1978) 

V. LEAD COMPUTER SCIENTIST 

Joel H. Saltz - Ph.D., Computer Science, Duke University, 1985. Parallel Computing with 
Emphasis on Systems and Algorithmic Issues. (Beginning July 1989) 

V. SENIOR STAFF SCIENTIST 

Gordon Erlebacher - Ph.D., Plasma Physics, Columbia University, 1983. Computational 
Fluid Dynamics. (November 1989 to November 1994) 
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Dimitri Mavriplis - Ph.D., Mechanical and Aerospace Engineering, Princeton University, 
1988. Grid Techniques for Computational Fluid Dynamics. (February 1987 to September 
1995) 

Piyush Mehrotra - Ph.D., Computer Science, University of Virginia, 1982. Assistant Pro- 
fessor, Department of Computer Science, Purdue University. Programming Languages for 
Multiprocessor Systems. (January to September 1992) 

Charles G. Speziale - Ph.D., Aerospace and Mechanical Sciences, Princeton University, 
1978. Fluid Dynamics with Emphasis on Turbulence Modeling and the Transition Process. 
(September 1987 to September 1992) 

Shlomo Ta’asan - Ph.D., Applied Mathematics, The Weizmann Institute of Science, 1985. 
Multigrid Methods for Partial Differential Equations. (July 1991 to July 1994) 

John R. Van Rosendale - Ph.D., Computer Science, University of Illinois, 1980. Parallel 
Systems and Algorithms. (July 1989 to July 1992) 

VI. SCIENTIFIC STAFF 

Fabio Bertolotti - Ph.D., Mechanical Engineering, Ohio State University, 1991. Stability 
Theory in Fluid Mechanics. (September 1991 to September 1993) 

Kurt M. Bryan - Ph.D., Mathematics, University of Washington, 1990. Theoretical and 
Computational Methods for Inverse Problems. (August 1990 to September 1993) 

Leon M. Clancy - B.S., Mechanical Engineering, University of Washington, 1971. System 
Manager. (December 1989 to Present) 

Thomas W. Crockett - B.S., Mathematics, College of William and Mary, 1977. Parallel 
Systems Research. (February 1987 to September 1992) 

Subhendu Das - M.S., Computer Science, College of William and Mary, 1990. Parallel Tools 
and Environments for Unstructured Scientific Computations. (June 1991 to June 1993) 

Thomas M. Eidson - Ph.D., Mechanical Engineering, University of Michigan, 1982. Parallel 
Techniques for Computational Fluid Dynamics. (August 1989 to August 1991) 

Yuh-Roung Ou - Ph.D., Aerospace Engineering, University of Southern California, 1988. 
Control Systems for Fluid Dynamics. (November 1988 to September 1991) 
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Peter W. Pretzel - Ph.D., Electrical Engineering, Technical University of Braunschweig, 
Germany, 1987. Reliability of Computing Systems. (March 1987 to September 1991) 

James Quirk - Ph.D., Computational Fluid Dynamics, Cranfield Institute of Technology, 
1991. Adaptive Methods for Partial Differential Equations. (June 1991 to June 1993) 

Sutanu Sarkar - Ph.D., Mechanical and Aerospace Engineering, Cornell University, 1988. 
Fluid Dynamics, Turbulence Modeling, Compressible Turbulence. (September 1988 to Septem- 
ber 1992) 

Sharon O. Scddougui - Ph.D., Applied Mathematics, University of Exeter, England, 1988. 
Compressible Fluid Dynamics. (June 1988 to August 1991) 

Ralph C. Smith - Ph.D., Numerical Analysis, Montana State University, 1990. Theoretical 
and Computational Issues Associated with Inverse Problems. (August 1990 to August 1993) 

Alan Sussman - Ph.D., Computer Science, Carnegie-Mellon University, 1991. Quantifying 
aixd Predicting Performance of Realistic Applications on Distributed Memory Parallel Com- 
puters. (August 1991 to August 1993) 

VII. VISITING SCIENTISTS 

Remi Abgrall - Ph.D., Mechanical Engineering, University of Paris, France, 1987. Research 
Scientist, INRIA (France). Computational Fluid Dynamics on Unstructured Meshes. (July 
1991 to August 1991) 

Shahid Bokhari - Ph.D., Electrical and Computer Engineering, University of Massachusetts, 
Amherst, 1978. Professor, Electrical Engineering, University of Engineering k Technology, 
Lahore, Pakistan. Parallel Computing Systems. (July 1991 to September 1991) 

Claudio Canuto - Ph.D., Mathematics, Turin, Italy, 1975. Researcher, Institute Analisi 
Numerica, Italy. Spectral Methods for Partial Differential Equations. (September 1991). 

Barbara Chapman - M.S., Mathematics, University of Canterbury, Christchurch, New Zealand, 
1985. Research Associate, Computer Science Department, University of Vienna. Compiler 
Development for Multiprocessors. (August 1991) 

William Criminale - Ph.D., Aeronautics, The Johns Hopkins University, 1960. Professor, 
Department of Applied Mathematics, University of Washington. Stability and Transition. 
(August 1991 to September 1991) 
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Kathleen E. Crowley - Ph.D., Computer Science, University of Washington, 1986. Sr. Pro- 
grammer, Department of Computer Science, Yale University. Environments for Parallel 
Computing Systems. (February 1991 to August 1991). 

Wai-Sun Don - Ph.D., Numerical Analysis, Brown University, 1989. Visiting Assistant Pro- 
fessor, Applied Mathematics, Brown University. Numerical Methods for Partial Differential 
Equations. (June 1991 to July 1991) 

Peter W. Duck - Ph.D., Fluid Mechanics, University of Southampton, United Kingdom, 
1975. Lecturer in Mathematics, Department of Mathematics, University of Manchester, 
United Kingdom. Numerical Solution of Unsteady Boundary Layer Equations. (July 1991 
to September 1991) 

Daniel J. Inman - Ph.D., Mechanical Engineering, Michigan State University, 1980. Associate 
Professor, Department of Mechanical and Aerospace Engineering, State University of New 
York at Buffalo. Numerical Methods for Inverse Problems in Distributed Systems. (June 
1991) 

Peter A. Jacobs - Ph.D., Computational Fluid Dynamics, University of Queensland, Aus- 
tralia, 1987. Compressible Fluid Dynamics. (January 1990 to October 1991) 

Mordechay Karpel - Ph.D., Aeronautics and Astronautics, Stanford University, 1980. As- 
sociate Professor, Department of Aerospace Engineering, Technion, Israel. Modeling of 
Aeroservoelastic Systems. (July 1991 to September 1991) 

FumioKojima- Ph.D., Control Theory, Kyoto University, Japan. 1985. Associate Professor, 
Department of Mechanical Engineering, Osaka Institute of Technology, Japan. Probabilistic 
and Stochastic Methods for Optimal Control Problems. (August 1991) 

Robert Krasny - Ph.D., Applied Mathematics, University of California, Berkeley, 1983. As- 
sociate Professor, Department of Mathematics, University of Michigan. Vortex Methods and 
Vortex Dynamics. (July 1991 - August 1991) 

Marten Landahl - Ph.D., Aeronautics, Royal Institute of Technology, Sweden, 1959. Profes- 
sor of Aeronautics and Astronautics, Massachusetts Institute of Technology. Transition and 
Turbulence. (August 1991). 

Jacques Liandrat - Ph.D., Fluid Mechanics, Uni versite Paul Sabatier, France, 1986. Research 
Scientist, IMST, France. Numerical Algorithms for Nonlinear PDE. (June 1991) 
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Moshe Matalon - Ph.D., Mechanical and Aerospace Engineering, Cornell University, 1978. 
Professor, Engineering Sciences and Applied Mathematics Department, Northwestern Uni- 
versity. Hydrodynamic Instability in Chemically Reacting Flows. (June 1991) 

K. W. Morton - Ph.D., Mathematics, New York University, 1964. Professor of Numerical 
Analysis, Computing Laboratory, Oxford University, England. Computational Fluid Dy- 
namics. (September 1991) 

H.S. Mukunda - Ph.D., Engineering Sciences, Indian Institute of Science, 1970. Professor, 
Department of Aerospace Engineering, Indian Institute of Science. High Speed Reacting 
Flows. (May 1991 to June 1991) 

Rolf Radespiel - Ph.D., Mechanical Engineering, Braunschweig Technical University, 1986. 
Head, Aerothermodynamics Section. Numerical Methods for High Speed Flows. (June 1991 
to September 1991) 

Asok Ray - Ph.D., Mechanical Engineering, Northeastern University, 1976. Professor, Me- 
chanical Engineering Department, The Pennsylvania State University. Control Systems for 
Flexible Space Structures. (May 1991 to June 1991) 

Philip L. Roe - Ph.D., Aeronautical Engineering, Cambridge, United Kingdom, 1962. Pro- 
fessor, Aerospace Engineering, University of Michigan, Ann Arbor. Computational Fluid 
Dynamics. (August 1991) 

Eitan Tadmor - Ph.D., Numerical Analysis, Tel-Aviv University, 1979. Professor, Depart- 
ment of Applied Mathematics, Tel-Aviv University. Numerical Methods for Partial Differ 
ential Equations. (July 1991 - August 1991) 

Siva Thangam - Ph.D., Mechanical Engineering, Rutgers University, 1980. Professor, De- 
partment of Mechanical Engineering, Stevens Institute of Technology. Computational Fluid 
Mechanics. (August 1990 to August 1991) 

Hein Tran - Ph.D., Mathematics, Renssela.er Polytechnic Institute, 1986. Assistant Professor, 
Department of Mathematics, North Carolina State University. Control Theory. (June 1991) 
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VIII. CONSULTANTS 


Loyce M. Adams - Ph.D., Applied Mathematics, University of Virginia, 1983. Associate Pro- 
fessor, Department of Applied Mathematics, University of Washington. Numerical Methods 
for Parallel Computing Systems. 

Dinshaw Balsara - Ph.D., Computational Fluid Dynamics, Astro Physics, University of 
Illinois at Urbana, 1990. Department of Physics and Astronomy, Johns Hopkins University. 
Parallel Implementation on Adaptive Godunov Schemes. 

Alvin Bayliss - Ph.D., Mathematics, New York University, 1975. Associate Professor, Tech- 
nological Institute, Northwestern University. Numerical Methods for Partial Differential 
Equations. 

Marsha J. Berger - Ph.D., Numerical Analysis, Stanford University, 1982. Research Asso- 
ciate, Courant Institute of Mathematical Sciences. Numerical Methods for Partial Differen- 
tial Equations. 

Harry Scott Berryman - B.S., Computer Science, Yale University, 1988. Graduate Fellow, 
Department of Computer Science, Yale University. Software for Parallel Systems. 

Percy Bobbitt - B.S., Aeronautics, Catholic University of America, 1949. NASA Langley 
Research Center - Retired. Fluid Mechanics. 

John A. Burns - Ph.D., Mathematics, University of Oklahoma, 1973. Professor, Virginia 
Polytechnic Institute and State University. Numerical Methods in Feedback Control and 
Parameter Estimation. 

Peter R. Eiseman - Ph.D., Mathematics, University of Illinois, 1970. Senior Research Sci- 
entist and Adjunct Professor, Department of Applied Physics and of Nuclear Engineering, 
Columbia University. Computational Fluid Dynamics. 

Robert E. Fennell - Ph.D., Mathematics, University of Iowa, 1969. Professor, Department 
of Mathematical Sciences, Clemson University. Control Theory for Multivariable Systems. 

Joel H. Ferziger - Ph.D., Nuclear Engineering, University of Michigan, 1962. Professor, 
Thermosciences Division, Department of Mechanical Engineering, Stanford University. Fluid 
Dynamics. 

Dennis B. Gannon - Ph.D., Computer Science, University of Illinois, 1980. Associate Pro- 
fessor, Department of Computer Science, Indiana University. Parallel Computation. 
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James F. Geer - Ph.D., Applied Mathematics, New York University, 1967. Professor, Sys- 
tems Science and Mathematical Sciences, Watson School of Engineering, Applied Science 
and Technology, SUNY-Binghamton. Perturbation Methods and Asymptotic Expansions of 
Solutions to Partial Differential Equations. 

Pratik Gupta - M.S., Computer Science, Virginia Commonwealth University, 1989. Systems 
Engineer, Simulation Associates Inc. /Taurus Technologies Inc. Parallel Computing Systems. 

Philip Hall - Ph.D., Mathematics, Imperial College, England, 1973. Professor, Department 
of Mathematics, University of Manchester, England. Computational Fluid Dynamics. 

Amiram Harten - Ph.D., Mathematics, New York University, 1974. Associate Professor, 
Department of Mathematics, Tel-Aviv University, Israel. Numerical Solution for Partial 
Differential Equations. 

Kazufumi Ito - Ph.D., Systems Science and Mathematics, Washington University, 1981. 
Assistant Professor, Department of Mathematics, University of Southern California. Control 
Theory. 

Thomas L. Jackson - Ph.D., Mathematics, Rensselaer Polytechnic Institute, 1985. Assistant 
Professor, Department of Mathematics, Old Dominion University. Numerical and Analytical 
Methods for Chemically Reacting Flows. 

Antony Jameson - Ph.D., Magnetohydro-Dynamics, Cambridge University, England, 1963. 
James S. McDonnell Distinguished Professor, Department of Mechanical and Aerospace 
Engineering, Princeton University. Computational Fluid Dynamics. 

Mark T. Jones - Ph.D., Computer Science, Duke University, 1990. Assistant Computer 
Scientist, MCS Division, Argonne National Labs. Parallel Algorithms for Numerical Linear 
Algebra. 

Harry F. Jordan - Ph.D., Physics, University of Illinois, 1967. Professor Department of Elec- 
trical and Computer Engineering, University of Colorado at Boulder. Parallel Computation. 

Ashwani K. Kapila - Ph.D., Theoretical and Applied Mechanics, Cornell University, 1975. 
Associate Professor, Department of Mathematical Sciences, Rensselaer Polytechnic Institute. 
Ordinary and Partial Differential Equations, Asymptotic Methods. 

David E. Keyes - Ph.D., Applied Mathematics, Harvard University, 1984. Assistant Pro- 
fessor, Mechanical Engineering, Yale University. Parallelizaton of Numerical Procedures 
Appropriate for the Study of Combustion. 
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Fumio Kojima - Ph.D., Control Theory, Kyoto University, Japan, 1985. Visiting Research 
Assistant Professor, Center for Applied Mathematical Sciences, University of Southren Cal- 
ifornia. Probabilistic and Stochastic Methods for Optimal Control Problems. 

Heinz-Otto Kreiss - Ph.D., Mathematics, Royal Institute of Technology, Sweden 1960. Pro- 
fessor, Department of Applied Mathematics, California Institute of Technology. Numerical 
Analysis. 

William D. Lakin - Ph.D., Applied Mathematics, University of Chicago, 1968. Professor, 
Department of Mathematics and Statistics, University of Vermont. Computational Fluid 
Dynamics. 

David G. Lasseigne - Ph.D., Applied Mathematics, Northwestern University, 1985. Assistant 
Professor, Department of Mathematics and Statistics, Old Dominion University. Computa- 
tional Fluid Dynamics. 

Anthony Leonard - Ph.D., Nuclear Engineering, Stanford University, 1963. Professor of 
Aeronautics, California Institute of Technology. Fluid Physics. 

Robert W. MacCormack - M.S., Mathematics, Stanford University. Professor, Department 
of Aeronautics and Astronautics, Stanford University. Computational Fluid Dynamics and 
Numerical Analysis. 

Ravi Mirchandaney - Ph.D., Computer Engineering, University of Massachusetts, 1987. Re- 
search Computer Scientist, Shell Oil Company, Houston, TX. Parallel Run-Time Support 
Systems. 

Seema Mirchandaney - M.S., Computer Science, University of Massachusetts, Amherst, 1990. 
Research Programmer, Science and Technology Center, Rice University. Parallel Program- 
ming Environments. 

Mark V. Morkovin - Ph.D., Applied Mathematics, University of Wisconsin, 1942. Profes- 
sor Emeritus, Department of Mechanical and Aerospace Engineering, Illinois Institute of 
Technology. Transition Process in Aerodynamics. 

Kirsten A. Morris - Ph.D., Electrical Engineering, University of Waterloo, 1989. Assistant 
Professor, Department of Applied Mathematics, University of Waterloo- Ontario, Canada. 
Control Theory. 

Nax)mi H. Naik - Ph.D., Mathematics, University of Wisconsin-Madison, 1987. Visiting 
Assistant Professor, Department of Mathematics, Vassar College. Multi-Grid Methods. 
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David M. Nicol - Ph.D., Computer Science, University of Virginia, 1985. Professor, Depart- 
ment of Computer Science, College of William and Mary. Mapping Algorithms onto Parallel 
Computing Systems. 

Roy A. Nicolaides - Ph.D., Mathematics, University of London, 1972. Professor, Depart- 
ment of Mathematics, Carnegie-Mellon University. Numerical Solution of Partial Differential 
Equations. 

James M. Ortega - Ph.D., Mathematics, Stanford University, 1962. Professor and Chairman, 
Department of Applied Mathematics, University of Virginia. Numerical Methods for Partial 
Differential Equations. 

Stanley J. Osher- Ph.D., Functional Analysis, New York University, 1966. Professor, Depart- 
ment of Mathematics, University of California at Los Angeles. Methods for the Numerical 
Analysis of Partial Differential Equations. 

Demetrius Papageorgiou - Ph.D., Mathematics, University of London, 1985. Assistant Pro- 
fessor, Department of Mathematics, New Jersey Institute of Technology. Computational 
Fluid Dynamics. 

Ugo Piomelli - Ph.D., Mechanical Engineering, Stanford University 1987. Professor, Depart- 
ment of Mechanical Engineering, University of Maryland. Subgrid Scale Reynold’s Stress 
Modelling and Large Eddy Simulation of Turbulent Flows. 

Terrence W. Pratt - Ph.D., Mathematics/Computer Science, University of Texas at Austin, 
1965. Professor, Department of Computer Science, University of Virginia. Programming 
Languages. 

Daniel A. Reed - Ph.D., Computer Science, Purdue University, 1983. Assistant Professor, 
Department of Computer Science, University of Illinois. Parallel Processing. 

Helen L. Reed - Ph.D., Engineering Mechanics, Virginia Polytechnic Institute and State 
University, 1981. Associate Professor, Department of Mechanical Engineering, Arizona State 
University. Computional Fluid Dynamics. 

Eli Reshotko - Ph.D., Aeronautics and Physics, California Institute of Technology, 1960. In- 
terim Dean, Case Western Reserve University, High Speed Aerodynamics with an Emphasis 
on Transition, Turbulence and Combustion. 
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Paul F. Reynolds - Ph.D., Computer Science, University of Texas at Austin, 1979. Assistant 
Professor, Department of Computer Science, The University of Virginia. Parallel Computing 
Systems. 

Paul E. Saylor - Ph.D., Mathematics, Rice University, 1986. Associate Professor, Computer 
Science Department, University of Illinois, Urbana. Iterative Solution of Linear Algebraic 
Equations and Algorithms for Parallel Computers. 

Jeffrey S. Scroggs - Ph.D., Computer Science, University of Illinois at Urbana, 1988. As- 
sistant Professor, Department of Mathematics, North Carolina State University. Domain 
Decomposition Techniques for Partial Differential Equations. 

Chi- Wang Shu - Ph.D., Mathematics, University of California, Los Angeles, 1986. Assistant 
Professor, Division of Applied Mathematics, Brown University. Partial Differential Equa- 
tions. 

Katepalli R. Sreenivasan - Ph.D., Aeronautical Engineering, Indian Institute of Science, 
1975. Professor and Chairman, Department of Mechanical Engineering, Yale University. 
Transition and Turbulence. 

Shlomo Ta’asan - Ph.D., Applied Mathematics, Weizmann Institute, 1984. Scientist, De- 
partment of Applied Mathematics, The Weizmann Institute of Science, ISRAEL. Multigrid 
Methods for Partial Differential Equations. 

Saleh Tanveer - Ph.D., Applied Mathematics, California Institute of Technology, 1984. Pro- 
fessor, Department of Mathematics, Ohio State University. Problems for Crystal Growth. 

Lu Ting - Ph.D., Aeronautics, New York University, 1951. Professor, Courant Institute of 
Mathematical Sciences, New York University. Fluid Mechanics. 

Eli Turkel - Ph.D., Applied Mathematics, New York University, 1970. Associate Profes- 
sor, Department of Applied Mathematics, Tel- Aviv University, Israel. Computational Fluid 
Dynamics. 

Bram van Leer - Ph.D., Theoretical Astrophysic, Leiden State University, The Netherlands, 
1970. Professor, Department of Aerospace Engineering, University of Michigan. Computa- 
tional Fluid Dynamics. 

Yun Wang - M.S., Applied Mathematics, University of Southern California, 1986. Research 
Associate, Center for Control Sciences, Division of Applied Mathematics, Brown University. 
Control Theory. 
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